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There  are  a wide  variety  of  electro-optical  (E-0)  sensors  employed 
in  the  modern  tactical  battlefield  environment.  Models  for  all  sources 
of  degradation  are  required  before  the  battlefield  performance  of  these 
E-0  sensors  can  be  predicted.  The  two  categories  of  degradation  sources 
are  the  ambient  (adverse  weather,  wind  driven  dust,  etc.)  and  the  battle- 
field induced  contaminants  (BIC,  any  adverse  sources  introduced  by  man). 
The  two  major  BIC  sources  are  countermeasures  smoke  and  the  dust  clouds 
due  to  tactical  munitions  impacting  the  ground.  In  this  report  we  begin 
the  task  of  modeling  the  munitions  dust  clouds. 

The  dust  models  ve  will  develop  will  be  the  first  generation  models. 
As  more  information,  primarily  experimental,  becomes  available,  succeed- 
ing improved  generations  of  models  can  be  developed.  Our  data  sources, 
both  theoretical  and  experimental,  are  severely  limited  at  this  time. 
Because  of  the  importance  of  tactical  dust  clouds,  a series  of  test 
programs  have  begun.  In  November  and  December  of  1977  six  dust  cloud 
tests  were  conducted  at  the  Dugway  Proving  Ground.  Dust  clouds  were 
generated  using  static  TNT  charges  to  simulate  the  detonation  of  high 
explosive  munitions  rounds.  Table  1 shows  trial  number,  the  munition 
simulated  and  the  amount  of  TNT  used.  Photographic  data  were  taken  of 
the  dust  clouds  produced.  Data  were  collected  from  various  electro- 
optical  instruments  operating  at  several  wavelengths.  At  Fort  Sill  in 
May  1978  tests  involving  live  firings  were  conducted.  Table  2 shows  the 
trials  and  the  number  of  rounds  fired  (both  single  and  multiple  round 
tests  were  conducted) . Again  both  photographic  and  electro-optical  data 
were  taken  of  the  resulting  dust  clouds.  In  October  1978  a series  of 
tests  of  both  static  and  live  rounds  were  conducted  at  White  Sands 
Missile  Range. 
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Trial  Number 


DPG-002-D1 

DPG-002-D2 

DPG-002-D4 

DPG-002-D3 

DPG-002-D5 

DPG-002-D6 


Simulated  Round 


81  mm  motar 

4.2  inci  mortar 

120  mm  mortar 

105  mm  tank  HEP 

155  mm  artillery  projectile 

8 inch  artillery  projectile 


Weight  of  TNT  (lbs) 


3.5 

5.0 

8.0 
6.75 

15.0 

35.0 


Table  2.  List  of  Fort  Sill  test  parameters 


Trial  Designation 

] 

Weapon  System 

T3 

155  mm,  M107 

T4-T8 

155  mm,  M107 

T9-T11 

155  mm,  M107 

T12-T20 

105  mm.  Ml 

T21 

105  mm,  Ml 

T22 

105  mm.  Ml 

Number  of  Rounds  Fired 


1 


Unfortunately  at  this  time  only  a limited  amount  of  reduced  data 
are  available  from  the  test  series  for  use  in  the  dust  cloud  modeling 
effort.  For  the  Dugway  tests,  a preliminary  copy  of  the  test  report 
was  made  available  to  us.  The  report  gives  the  data  taken  with  the 
electro-optical  sensors.  In  addition  the  dust  cloud  length,  width, 
and  vertical  extent  at  10  second  intervals  (out  to  about  60-90  seconds) 
were  separately  provided.  For  the  Fort  Sill  tests,  some  preliminary 
data  and  the  cloud  dimensions  (out  to  about  15  seconds)  were  provided. 
No  data  were  available  from  the  just  completed  White  Sands  Missile 
Range  tests. 
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Figure  1 is  a simplified  flow  chart  of  the  elements  required  to 
construct  a model  for  the  effects  of  dust  clouds  in  a tactical  battle- 
field environment.  The  primary  purpose  of  this  report  is  to  develop 
the  models  for  the  dust  cloud  itself  (those  items  in  the  CLOUO  MODEL 
box  of  Figure  1).  We  will  develop  models  for  the  munition  dust  cloud 
size,  location  and  distribution  of  moss  within  the  cloud  as  a function 
of  time.  For  the  cloud  model  development,  it  was  necessary  to  also 
develop  the  MIL  MODIiL  parameters  and  parts  of  the  PARTICLE  MODEL, 
RADIATION  TRANSPORT  and  OUTPUTS  parameters.  The  basic  modeling  philoso- 
phy is  to  use  first  principle  physics  or  similarity  theory  to  set  the 
form  of  the  equation  for  a given  effect,  then  use  the  experimental  data, 
if  they  exist,  to  set  the  magnitudes  of  any  adjustable  parameters  in  the 
equation.  We  first  give  an  overview  of  the  whole  dust  cloud  model.  We 
then  develop  the  details  of  the  model  in  the  following  sections.  The 
MIE  MODEL  details  are  given  in  the  appendix. 

Figure  2 is  a sketcli  or  the  dust  cloud  formation  and  development. 

We  ignore  the  extremely  fast  formation  processes  and  assume  an  initial 
cloud  is  formed  instantaneously  at  burst  time.  We  assume  that  the  dust 
mass  lofted  by  the  munition  explosion  consists  of  dust  particles  ranging 
in  size  from  0 to  1 cm  in  diameter.  The  upper  limit  of  particle  sizes 
is  somewhat  arbitrary  but  is  of  little  significance  since  the  larger 
particles  fall  out  so  rapidly. 

We  divide  the  particle  size  range  into  a number  of  individual  size 
groups.  Table  3 shows  a preliminary  choice  of  size  groups.  We  want  a 
large  enough  number  of  size  groups  to  provide  a continuous  cloud  descrip- 
tion, but  a small  enough  number  for  rapid  computer  calculations.  All 
size  groups- are  assumed  to  occupy  the  same  volume  of  the  initial  cloud 
at  time  zero.  We  model  the  rise,  expansion,  wind  transport  and  diffu- 
sion of  each  size  group  separately.  The  lightest,  particles  will  rise, 
expand,  transport  and  diffuse  at  the  fastest  rates,  while  the  heavier 
particles  lag  behind  and  fall  out.  Figure  3 shows  a sketch  of  the  size 
groups  in  the  dust  cloud. 
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Figure  2.  Munition  dust  cloud  development. 
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Table  3.  Particle  size  groups 


Group  Range  of  Particle  Diameters  (cm) 


0 - 0.001 
0.001  - 0.002 
.002  - 
0.004  - 0.006 
0.006  - 0.008 
0.008  - 0.01 
0.01  - 0.02 
.02  - 0.03 
0.03  - 0.04 
0.04  - 0.05 
0.05  - 0.06 
0.06  - 0.07 
0.07  - 0.08 
0.08  - 0.09 
0.09  - 0.1 


- 0.4 

- 0.6 
- 0.8 
- 1.0 


Consider  an  x,  y,  z rectangular  coordinate  system  with  origin  at 
ground  level  at  the  explosion  point,  z is  in  the  vertical  direction; 
x is  along  the  wind  direction;  and  y is  in  the  cross  wind  direction. 

We  take  the  geometry  of  each  size  group  to  be  -spheroidal  with  cloud 

radii  R , R and  R (see  Figure  4).  The  radii  are  measured  from  the 
x y z 

centroid  of  each  size  group,  (x  , yQ,  zc).  We  will  model  the  centroid 
positions  and  radii  as  a function  of  time. 

Each  size  group  has  a relatively  small  range  of  particle  sizes  in  it. 
We  assume  no  size  separation  within  a given  group.  That  is,  each  element 
of  mass  anywhere  in  the  volume  of  a given  size  group  has  the  same  particle 


SMALLEST  SIZE  GROUP 


/ 


v, 


LARGEST  SIZE 


Figure  3.  Size  groups  in  dust  cloud. 


Figure  4.  Spheroidal  geometry. 
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size  distribution.  This  uniform  size  mixing  assumption  allows  us  to 

calculate  the  scattering  and  absorption  mass  coefficients  for  each  size 

group  only  once.  The  scatter  and  attenuation  of  each  size  group  is 

then  only  a function  of  the  mass  penetrated.  We  assume  that  within  each 

size  group  spheroid,  the  mass  distribution  is  Gaussian. 

2 2 2' 
r liY-x  i iV-Vv  / Z-Z 

M 


P(x,y,z)  = 


3/2 

(2*>5/  Vy°z 


exp 


h|R)-R)-R)|] 


(1) 


where 


= dust  mass  density  at  point  (x,y,z)  within  me  size  group 
volume 


(R  m"3) 


M 


total  mass  in  given  size  group  (g) 


2 2 
Ox  = variance  of  mass  density  in  the  x direction  (m  ) 


x£  <=  x-coordinate  of  centroid  of  size  group  volume 


Cm) 


M is  a constant  determined  by  the  dust  loading  model.  The  variances 
and  the  centroid  coordinates  are  functions  of  time.  We  assume  the 
visible  cloud  radii  and  the  mass  variances . arc  related  by 


R.  = C O. 
l r i 


(2) 


where 


l - x,  y or  z 


Cr  = constant 


The  value  of  Cr  depends  upon  the  measurement  technique.  The  visible 


cloud  dimensions  are  generally  determined  by  photographic  data.  If  we 
assume  the  image  brightness  at  any  point  is  pi'oportional  to  the  mass 
penetrated  through  the  cloud  at  that  point,  then  a common  assumption  is 
the  visible  cloud  edge  occurs  at  a mass  penetrated  value  which  is  a 
constant  fraction  of  that  at  the  cloud  center.  A typical  fraction  value 
is  0.1,  in  which  case  Cr  is  2.15  (Reference  5).  Other  assumptions  and 
measurement  techniques  give  different  values  for  C^,  but  is  almost 
always  in  the  range  from  1 to  3.  Until  we  have  more  data  for  the  muni- 
tion clouds,  we  take 


Cr  = 2.15 


(3) 
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SECTION  2 

INITIAL  DUST  CLOUD 
RADIUS  AND  SHAPE 


: 


Consider  a munitions  round  impacting  the  ground.  Depending  on  the 
fusing,  the  round  can  detonate  at  or  under  the  ground  surface.  The  ex- 
plosive energy  release  will  lead  to  crater  formation  and  the  lofting  of 
dust  into  the  air.  Part  of  the  energy  released  will  couple  into  the  air 
and  appear  as  hydrodynamic,  ie,  as  shock,  energy.  The  formation  of  a 
blast  shock  will  lead  to  a very  rapid  expansion  of  the  explosive  cloud. 

At  first  the  pressure  inside  the  cloud  is  much  greater  than  the  ambient 
pressure.  But  as  the  cloud  expands,  the  internal  pressure  drops.  When 
the  cloud  comes  into  approximate  pressure  equilibrium  with  the  ambient 
atmosphere,  the  rapid  expansion  ceases  and  the  initially  high  internal 
air  flow  velocities  rapidly  decrease. 

The  initial  size  of  the  dust  cloud  is  not  a precisely  defined  quantity. 
The  early  rapid  expansion  slows  and  continuously  merges  into  the  slow 
expansion  of  the  rise  and  diffusion  phase.  The  time  scale  of  the  air 
shock  motion  is  extremely  short;  pressure  equilibrium  is  reached  in  a 
few  milliseconds  or  less.  However,  it  may  take  several  tenths  of  sec- 
onds for  the  initial  high  dust  particle  velocities  to  slow.  We  take  the 
initial  cloud  size  as  the  size  at  the  end  of  the  rapid  expansion  phase 
and  the  beginning  of  the  slow  expansion  phase.  We  ignore  the  details  of 
the  rapid  expansion  phase  and  assume  the  initial  dust  cloud  is  formed 
instantaneously  at  burst  time. 

We  will  rely  upon  the  theory  of  shock  physics  for  the  functional  form 
of  the  initial  dust  cloud  radius.  We  then  use  the  experimental  data  to 
Set  the  value  of  any  adjustable  parameters.  The  motion  of  a spherical 
shock  in  uniform  air  with  a constant  gas  gamma  can  be  expressed  in  similarity 
equations.  The  normalizing  distance  and  time  constants  are 


1C 


Kjj  = hydrodynamic  yield  of  the  explosion  producing 
the  shock  (ergs) 

-2 

P «=  ambient  pressure  (dynes  cm  ) 

_3 

p «=  ambient  density  (g  cm  ). 

Then  in  terms  of  R/R  , t/tQ  all  shocks  are  similar.  Figure  5 from 
Reference  1 gives  the  properties  at  the  shock  front.. 

We  assume  that  the  initial  radius  of  the  cloud  is  the  shock  radius 
at  pressure  equilibrium  time.  From  hydro  code  calculations  and  observa- 
tions of  nuclear  fireballs,  the  pressure  equilibrium  time  is  approximately 


0.162 


at  which  time 


0.562 


If  we  assume  1 Lb  TNT  = S x 105  cal  = 2.00  x 1013  ergs,  and 


"p~  s (by  custom  in  weapon  effects  scaling,  p/pQ  is 

° 0 used  instead  of  P/PD) 

Then  the  initial  explosive  cloud  radius  can  be  written 


For  a free  air  explosion,  f would  be  the  fraction  of  the  total  yield 
which  appears  as  effective  hydrodynamic  energy.  For  a munition  explosion 
in  free  air, 

f « 1 

For  a conventional  nuclear  device  detonated  in  free  air, 

f = 0.633 

The  remaining  fraction  of  the  nuclear  yield  is  radiated  away  from  the 
nuclear  fireball  as  thermal  emissions  (an  effect  not  present  in  the 
munitions  explosion). 

For  our  case  of  surface  or  below  surface  munitions  detonations,  f 
can  be  considered  an  empirical  factor,  to  be  set  primarily  from  the 
experimental  data,  f will  depend  on  shell  type,  explosion  depth,  and 
soil  type.  Table  4 from  Reference  2 gives  the  currently  available 
dust  cloud  size  data  from  the  Dugway  dust  tests.  The  data  arc  given  at 
10-second  intervals,  with  initial  dimensions  given  at  z-time  (which  we 
assume  to  be  a time  very  shortly  after  detonation).  The  initial  Dugway 
clouds  arc  roughly  spherical.  We  take  as  the  experimental  initial  radius 
one  half  (the  test  data  are  diameters)  of  the  average  of  the  X,  Y,  and 
Z initial  dimensions.  Figure  6 shows  the  initial  radius  versus  yield 
and  the  fit  to  the  data.  The  Dugway  explosives  were  static  TNT  charges 
buried  so  that  the  top  of  the  charge  was  at  the  ground  surface.  For 
these  type  of  charges  and  for  the  Dugway  soil  type, 

f = 0.84  (9) 

Tables"  5 through  13  from  Reference  3 give  the  currently  available 
data  for  single  round  tests  from  the  Fort  Sill  test  series.  Trial  T-3 
is  a 155-mm  round;  the  other  eight  tests  are  for  the  105-nun  round.  The 
size  data  are  given  at  0.1  second  intei'vals  for  the  first  second,  then 
at  1 second  intervals  thereafter.  There  is  considerable  scatter  in  the 
size  data  from  one  test  to  the  next.  But  in  general  the  early  cloud 
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Table  4.  Dugway  dust  cloud  measurements 


All  data  in 

meters;  X - 

X Y 

length, 

Z X 

Y = width, 

Y Z 

Z = height 

X Y 

Z 

X 

Y 

Z 

Z Time 

01 

3 5 

3 

5 

D2 

8 3 

D3 

6 6 

6 

D4 

4 4 

6 

Z + 10  sec. 

24  -22 

3 

40 

22  15 

20  19 

7 

39 

20 

14 

Z + 20  " 

58  -28 

3 

64 

23  25 

87  38 

11 

67 

34 

18 

Z + 30  " 

50  23 

13 

66 

25  34 

99  46 

21 

88 

38 

27 

Z + 40  " 

51  20 

13 

57 

40  30 

143  50 

24 

84 

53 

25 

Z + 50  " 

47  19 

13 

65 

39  31 

Cloud 

78 

58 

23 

Z + 60  " 

46  18 

13 

41 

48  30 

moved 

100 

64 

28 

Z + 70  " 

Z + 80  " 

Z + 90  " 

Z + 100  " 

z + no  * 

7 + " 

No  visible 
cloud 

No  visible 
cloud 

outside 

camera 

coverage 

Cloud 

moved 

outside 

camera 

coyerage 

X 

Y 

Z 

X 

Y 

Z 

D5 

06 

Z Time 

10 

N* 

9 

12 

N* 

11 

Z * 10 

sec. 

42 

0 

12 

56 

0 

26 

Z + 20 

m 

50 

12 

71 

D 

38 

Z + 30 

•i 

47 

0 

19 

69 

A 

53 

Z + 40 

it 

58 

A 

10 

66 

T 

49 

Z + 50 

n 

60 

T 

17 

84 

A 

47 

Z + 60 

H 

51 

A 

18 

ClOUd 

moved 

Z + 70 

H 

53 

12 

outside 

camera 

Z + 80 

H 

55 

16 

coverage 

Z + 90 

it 

59 

15 

*Dust  from  high  winds  obscured  view  of  camera.  Width  measurements 
cannot  be  provided. 
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Table  5.  Initial  cloud  dimensions  through  time  for  Trial  T-3. 


TIME 

LENGTH 

(meters) 

WIDTH 

(meters) 

HEIGHT 

(meters) 

954:01 

10 

ND 

3 

954:01.1 

13 

ND 

3 

954:01.2 

14 

ND 

3 

954:01.3 

15 

ND 

3 

954:01.4 

17 

ND 

4 

954:01.5 

17 

ND 

4 

954:01.6 

18 

ND 

4 

95^,4:01.7 

19 

ND 

4 

954:01.8 

19 

ND 

4 

954:01.9 

19 

ND 

4 

954:02 

20 

ND 

5 

954:03 

22 

ND 

5 

954:04 

26 

ND 

6 

954:05 

30 

ND 

9 

954:06 

30 

ND 

10 

954:07 

30 

ND 

10 

954:08 

33 

ND 

11 

954 : 09 

35 

ND 

12 

954:10 

43 

ND 

13 

954:11 

48 

ND 

13 

ND:  Width  of  cloud  not  available.  Cloud  doesn't  come  into  camera’s 


field  of  view  until  954:13. 


Table  6.  Initial  cloud  dimens 


TIME 

LENGTH 

(meters) 

WIDTH 

(meters) 

HEIGHT 

(meters) 

921:00 

4 

4 

2 

921:00.1 

4 

9 

3 

921:00.2 

7 

9 

3 

921:00.3 

7 

11 

3 

921:00.4 

8 

11 

3 

921:00.5 

8 

11 

4 

921:00.6 

8 

11 

4 

921:00.7 

9 

12 

4 

921:00.8 

10 

12 

5 

921:00.9 

10 

12 

6 

921:01 

11 

12 

6 

921:02 

15 

16 

7 

921:03 

17 

18 

8 

921:04 

17 

22 

9 

921:05 

18 

23 

9 

921:06 

18 

25 

9 

921:07 

19 

25 

10 

92,1:08 

21 

33 

12 

921:09 

21 

36 

12 

921:10 

23 

37 

12 

921:11 

26 

39 

11 

921:12 

36 

ND 

11 

921:13 

42 

ND  ! 

11 

921:14 

42 

ND 

10 

921:15 

42 

ND 

10 

ND:  Width  of  cloud  not  available.  Camera  stopped  at  921:11, 
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HEIGHT 

(meters) 


LENGTH 
(meters ) 


TIME 


927:02  1 
927:02.1  3 
927:02.2  5 


927:02.3  7 
927:02.4  7 
927:02.5  8 
927:02.6  8 
927:02.7  8 

92V: 02. 8 8 


927:02.9  9 
927:03  9 
927:04  11 
927:05  13 
927:06  15 
927:07  19 
927:08  21 
927:09  26 


Table  7.  Initial  cloud  dimensions  through  time  for  Trial  T-13. 


Table  9.  Initial  cloud  dimensions  through  time  for  Trial  T-15. 


957:59 

957:59.1 

957:59.2 

957:59.3 

957:59.4 

957:59.5 

957:59.6 

957:59.7 

95^7: 59.8 

957:59.9 

958:00 

958:01 

958:02 

958:03 

958:04 

958:05 

958:06 

958:07 

958':  08 

958:09 

958:10 

958:11 

958:12 

958:13 

958:14 


WIDTH 

(meters) 


HEIGHT 

(meters) 


Table  10.  Initial  cloud  dimensions  through  time  tor  Trial  T-16. 


TIME 

LENGTH 

(meters) 

WIDTH 

(meters) 

HEIGI 

(metei 

1004:01 

ND 

2 

3 

1004:01.1 

NO 

3 

3 

1004:01.2 

ND 

3 

3 

1004:01.3 

ND 

5 

3 

1004:01.4 

ND 

5 

3 

1004:01.5 

ND 

6 

3 

1004:01.6 

ND 

7 

3 

1004:01.7 

ND 

7 

3 

10t)4:01.8 

ND 

7 

3 

1004:01.9 

ND 

7 

3 

1004:02 

ND 

7 

4 

1004:03 

ND 

8 

5 

1004:04 

ND 

10 

6 

1004:05 

ND 

13 

7 

1004:06 

ND 

18 

7 

1004:07 

ND 

19 

7 

1004:08 

ND 

23 

7 

1004:09 

ND 

23 

7 

1004: 10 

ND 

23 

7 

1004:11 

ND 

23 

7 

1004:12 

ND 

28 

8 

1004:13 

ND 

28 

9 

1004:14 

ND 

28 

10 

1004:15 

ND 

29 

10 

1004:16 

ND 

31 

11 

1004:17 

ND 

31 

11 

1004:18 

ND 

31 

13 

ND:  Length  of  cloud  not  available.  Camera  did  not  run. 
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Table  11.  Initial  cloud  dimensions  through  time  for  Trial  T-17 


1010:01 

1010:01.1 

1010:01.2 

1010:01.3 

1010:01.4 

1010:01.5 

1010:01.6 

1010:01.7 

1010:01.8 

1010:01.9 

1010:02 

1010:03 

1010:04 

1010:05 

1010:06 

1010:07 

1010:08 

1010:09 

1010:10 

1010:11 

1010:12 

1010:13 

1010:14 

1010:15 

1010:16 

1010:17 


LENGTH 

(meters) 


HEIGHT 

(meters) 


Length  of  cloud  not  available.  Camera  did  not  run 


msm 
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Table  12.  Initial  cloud  dimensions  through  time  for  Trial  T-19. 


TIME 

WIDTH 

(meters) 

HEIGHT 

(meters) 

1141:00 

ND 

5 

1 

1141:00.1 

ND 

7 

2 

1141:00.2 

NO 

8 

2 

1141:00.3 

NO 

8 

2 

1141:00.4 

ND 

9 

2 

1141:00.5 

ND 

9 

3 

1141:00.6 

ND 

10 

3 

1141:00.7 

ND 

10 

3 

1141:00.8 

ND 

10 

3 

1141:00.9 

ND 

10 

3 

1141:01 

ND 

10 

3 

1141:02 

ND 

13 

3 

1141:03 

ND 

15 

3 

1141:04 

ND 

15 

5 

1141:05 

ND 

16 

6 

1141:06 

ND 

18 

6 

1141:07 

ND 

19 

7 

'-'1141:08 

ND 

19 

7 

1141:09 

ND 

23 

8 

1141:10 

ND 

23 

8 

ND:  Length  of  cloud  not  available.  Camera  did  not  run 
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Table  13.  Initial  cloud  dimensions  through  time  for  Trial  T-20 


TIME 

LENGTH 

(meters) 

WIDTH 

(meters) 

HEIGHT 

(meters) 

1148:00 

ND 

6 

2 

• 

1148:00.1 

NO 

6 

2 

• 

1148:00.2 

ND 

6 

2 

1148:00.3 

ND 

8 

2 

1148:00.4 

NO 

10 

2 

1148:00.5 

NO 

11 

2 

1148:00.6 

ND 

11 

2 

1148:00.7 

ND 

12 

3 

1148:00.8 

ND 

12 

3 

1148:00.9 

ND 

12 

3 

1148:01 

NO 

15 

3 

1148:02 

b 

15 

3 

1148:03 

ND 

15 

5 

1148:04 

ND 

16 

5 

1148:05 

ND 

18 

6 

1148:06 

NO 

19 

8 

NO:  Length  of  cloud  not  available.  Camera  did  not  run. 


30 


. 


r 

& 

i"  ■ 


. 


(time  less  than  about  0.3  seconds)  is  not  spherical;  the  vertical  dimension 
is  about  half  of  the  horizontal  dimension.  In  contrast  to  the  static 
charges  of  the  Dugway  tests,  the  Fort  Sill  tests  are  firings  of  live 
rounds.  The  rounds  impact  the  ground  at  an  angle  with  considerable 
velocity,  penetrate  the  ground  to  the  fusing  depth,  and  explode.  In  a 
dynamic  detonation  of  this  sort,  it  is  not  surprising  that  the  result- 
ing dust  cloud  is  nonspherical . The  forward  momentum  of  the  shell  may 
lead  to  an  initial  cloud  that  is  elongated  along  the  shell  track  with  a 
suppressed  vertical  dimension.  At  present  the  data  are  too  sketchy  to 
draw  any  firm  conclusions.  We  will  assume  that  the  initial  cloud  radius 
in  the  horizontal  direction  perpendicular  to  the  shell  track  is  given  by 
the  pressure  equilibrium  radius.  Equation  8.  The  initial  dimensions 
in  the  horizontal  along  track  direction  and  the  vertical  direction  will 
be  taken  as  constants  times  the  cross  track  dimension. 

The  preliminary  Fort  Sill  data  give  length,  width,  and  height  meas- 
urements for  the  dust  cloud,  but  unfortunately  give  neither  the  shell 
track  direction  nor  the  wind  direction.  At  early  times  we  expect  the 
cloud  to  be  elongated  along  the  shell  track  direction  and  at  later 
times  elongated  along  the  wind  direction.  We  use  the  width  dimensions 
of  Tables  9 through  13  for  both  the  initial  cross  track  cloud  size  and 
the  later  cross  wind  dimension.  The  cross  track  dimension  expands 
fairly  rapidly  for  the  first  0.2  to  0.3  seconds,  then  expands  at  a 
considerably  slower  rate  thereafter.  We  take  as  the  initial  size  of 
the  dust  cloud  the  width  radius  at  about  0.3  seconds,  corrected  for  the 
amount  of  expansion  due  simply  to  cloud  rise  in  this  time.  For  a 
munition  explosion,  we  expect  f to  have  its  largest  value,  f = 1,  for 
a free  air  burst  and  then  decrease  as  the  depth  of  burial  of  the  burst 
increases.  To  determine  f for  the  Fort  Sill  munitions,  we  utilize 
both  the  initial  size  data  and  the  later  time  expansion  data  (detailed 
in  the  next  section).  A reasonable  fit  to  both  early  and  late  data 
gives  for  the  105-  and  155-mm  rounds  and  the  Fort  Sill  soil  type 

f = 0.5  (10) 
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We  assume  the  initial  shape  of  any  dust  cloud  is  spheroidal,  with 
initial  radii 


Rj,  = initial  horizontal  radius  along  the  shell  track  (m) 
Rp  = initial  horizontal  radius  perpendicular  to  the  shell 
track  (m) 

Ry  = initial  vertical  radius  (m) . 

We  take 


CT  RI 


RP  = CP  RI 


Rv  = Ri 


(lla) 

(lib) 

(11c) 


where  Rj  is  the  initial  radius  given  by  Equation  8,  and  C.p,  Cp,  and 
Cy  are  shape  constants.  The  shape  constants  depend  upon  munition  and 
soil  type.  For  the  static  charges  and  the  soil  type  at  Dugway  Proving 
Grounds 


= ^P  = Cy  = ^ 


(12) 


For  the  105-  and  155-mm  rounds  and  the  soil  types  at  Fort  Sill 


CT«  2 


Cp  = 1 


C = — 

Lv  2 


(13a) 

(13b) 

(13c) 


We  have  reasonable  confidence  in  the  Cp  and  Cy  values  for  the  Fort 
Sill  shape  factors;  see  the  early  time  fits  of  Figures  10  and  11 


The  CT  value  is  less  certain.  The  chosen  value  of  2 is  not  incompatible 


with  the  width  to  length  ratios  of  Tables  6 to  8 for  the  first  several 
tenths  of  seconds.  Moreover  the  value  2 conserves  the  total  volume; 
that  is,  the  spheroid  volume  is  the  same  as  the  volume  of  a sphere 
with  the  shape  factors  equal  to  unity.  As  yet  we  do  not  have  sufficient 
data  to  take  the  conservation  of  volume  as  a general  rule. 


- - T . 
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As  more  data  become  available  for  different  munitions  and  soils, 
general  rules  can  be  formulated  for  the  functional  dependencies  of  f 
and  the  shape  factors.  The  initial  cloud  size  and  shape  are  significant 
only  for  the  first  several  seconds.  The  rapid  expansion  rates  in  the 
wind,  cross  wind,  and  vertical  directions  quickly  dominate  any  initial 
effects.  The  expansion  rates  determine  the  cloud  size  and  shape  and 
any  initial  influence  is  quickly  lost. 


SECTION  3 

DUST  CLOUD  RISE  PHASE 


The  initial  explosively  formed  cloud  will  have  an  average  internal 
temperature  higher  than  the  ambient  temperature.  Consequently,  the 
cloud  will  rise  buoyantly.  The  early  cloud  is  highly  turbulent.  The 
differential  velocity  between  the  rising  cloud  and  the  ambient  air  is 
a continuing  source  of  cloud  turbulence  during  the  rise  phase.  The 
initial  and  continuing  turbulence  during  rise  will  cause  the  cooler 
ambient  air  to  be  entrained  and  mixed  with  the  warmer  cloud  air.  The 
cloud  will  continue  rising  until  the  cooling  due  to  mass  entrainment 
and  adiabatic  expansion  bring  the  cloud  into  temperature  equilibrium 
with  the  atmosphere.  In  this  section  we  will  model  the  cloud  size, 
shape,  and  altitude  time  history  during  the  rise  phase. 

For  experimental  data  we  have  the  cloud  dimensions  given  in  Tables 
4 through  15  from  the  Dugway  and  Fort  Sill  test  series.  Unfortunately 
at  this  early  stage  in  the  test  data  reduction,  there  are  no  cloud 
altitude  data  available.  We  will  proceed  in  the  modeling  of  the  dust 
cloud  rise  phase  as  follows.  We  will  use  shock  physics  to  determine 
the  temperature  of  the  initial  cloud.  We  then  use  a theoretical  buoyant 
bubble  rise  with  entrainment  model  to  calculate  the  cloud  size  and 
altitude  time  history.  We  then  set  any  adjustable  parameters  to  give  a 
good  match  to  the  experimental  size  data. 

The  shock  code  of  Reference  1 was  used  to  calculate  the  tempera- 
ture, density,  and  pressure  profiles  inside  a spherical  shock  in  free 
air  as  a function  of  time.  At  pressure  equilibrium  time,  the  shock 
front  is  at  the  pressure  equilibrium  radius  of  liquation  8 . The 
average  temperature  as  defined  by 
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T = 1.4  T °K  (15) 

o 

where  Tq  is  the  ambient  temperature.  As  time  increases,  the  shock 
front  moves  farther  out;  the  density  inside  the  Rj  radius  decreases 
to  an  average  value 


but  the  average  temperature  remains  very  close  to  1.4  T . Hence  we 
take  the  temperature  of  the  initial  dust  cloud  to  be  1.4  TQ;  this  cor- 
responds to  403°K  for  a typical  ambient  temperature  of  288°K. 

We  will  use  the  GE — TEMPO  buoyant  fireball  rise  computer  model  to 
calculate  the  rise  and  expansion  of  the  dust  cloud.  Details  of  the 
rise  model  are  given  in  Reference  4 ; for  reader  convenience  we  give 
a very  brief  synopsis.  The  basic  rise  model  assumptions  are: 

• The  cloud  is  in  pressure  equilibrium  with  the  ambient 
atmosphere  at  all  times  during  the  rise  history. 

• The  cloud  is  isothermal  at  all  times. 

• As  the  cloud  rises,  it  entrains  a constant  fraction,  a,  of 
the  ambient  air  encountered.  The  entrained  air  is  instantly 
mixed  with  the  cloud  air. 

• The  drag  force  on  the  rising  cloud  can  be  described  by  a 
constant  drag  parameter,  Cp. 

• The  shape  of  the  rising  cloud  does  not  change;  the  cloud 
merely  gets  larger. 

With  these  assumptions  the  equations  for  the  cloud  temperature  rise 
velocity,  and  total  mass  are 


i 
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where 


temperature  of  cloud  (°K) 

_2 

pressure  of  cloud  (dynes  cm  ) 
cloud  rise  velocity  (cm  s~*) 

2 

cross  sectional  area  (perpendicular  to  v)  (cm  ) 
mass  of  air  in  cloud  (gm) 

number  of  moles  of  hot  air  of  temperature  T and  pressure 
P per  mole  of  cold  air 

effective  gas  constant  at  temperature  T and  pressure  P 
constant  entrainment  parameter 
constant  drag  parameter 

_2 

acceleration  of  gravity  (cm  s ) 

’ , Z , Y_  are  the  corresponding  quantities  for  the 
ambient  air. 


In  Equation  16  the  first  term  on  the  right-hand  side  is  the 
adiabatic  expansion  term,  the  second  term  is  the  mass  entrainment  term, 
and  the  third  and  fourth  terms  are  the  changing  equation  of  state  of 
air  terms.  Tho  third  and  fourth  terms  are  significant  for  the  high 
temperatures  of  a nuclear  fireball  cloud,  but  are  insignificant  for  the 


R 
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lower  temperatures  of  a munitions  cloud.  The  shape  of  the  cloud  is 
specified  by  the  user — whether  spherical,  spheroidal,  toroidal,  or 
other.  Given  the  initial  cloud  size  and  temperature,  the  rise  program 
numerically  integrates  the  rise  equations  forward  in  time  and  prints 
out  the  cloud  altitude,  size,  and  temperature. 

The  important  parameters  controlling  the  rise  of  a buoyant  cloud 
are  the  hydrodynamic  yield,  flV,  which  determines  the  initial  size  and 
temperature  of  the  cloud,  -md  the  mass  entrainment  factor,  a.  The  drag 
coefficient,  C^,  has  no  effect  upon  the  cloud  size  and  temperature 
versus  altitude,  but  merely  affects  the  time  at  which  the  cloud  reaches 
that  altitude.  The  present  rise  and  expansion  experimental  data  base 
is  incomplete;  so  the  cloud  shape,  entrainment  factor,  and  drag  coef- 
ficient cannot  be  determined  with  a high  confidence  level.  We  again 
employ  the  reasonable  procedure  of  using  the  theoretical  rise  program 
to  determine  the  functional  form  of  the  cloud  rise  and  expansion;  we 
then  adjust  the  resulting  equations  to  fit  the  existing  experimental 
data  (in  this  case,  the  expansion  data  since  there  are  no  rise  data). 

We  choose  the  simplest  case  for  the  rise  program  calculation,  a spherical 
cloud.  We  take  hydrodynamic  yields  (fW)  of  1,  10,  100,  and  1000  pounds 
of  TNT,  and  a range  of  entrainment  factors  (a)  of  0.5,  1,  and  2.  The 
drag  coefficient  (C^)  is  set  equal  to  0.8  (a  typical  value  for  low- 
altitude  nuclear  clouds).  The  initial  cloud  radius  is  calculated 
from  Equation  8 and  the  initial  cloud  temperature  was  set  at  403°K 
(1.4  times  the  ambient  temperature  of  288°K).  Figure  7 shows  the 
altitude  of  the  cloud  center  as  a function  of  time  for  a neutral 
stability  atmosphere  (see  Post  Rise  Phase  section).  Figure  8 
shows  the  cloud  radius  versus  time.  In  the  two  figures,  the  point  at 
which  the  sqlid  and  the  dotted  lines  join  is  the  point  at  which  the 
cloud  temperature  has  fallen  to  about  5°K  above  the  ambient  temperature. 

The  idealized  turbulent  entrainment  of  ambient  air  by  the  rising 
spherical  cloud  results  in  a radial  growth  which  is  described  by  the 
classical  diffusion  equation 
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Ks(t) 


= \//2)Tt~rKr 

s so 


m 


(19) 


where 


l<s(t)  = radius  of  the  rising  spherical  cloud  at  time  t (m) 


2 -1, 


effective  turbulent  radial  diffusion  coefficient  (m  s ) 


Rso  = initial  radius  of  sphere  (m) 

Rso  is  calculated  from  the  pressure  equilibrium  radius  Equation  8. 

An  excellent  fit  for  K is 

s 


K = 
s 


a >/fW 


2 -1 
m s 


(20) 


The  altitude  of  the  center  of  the  rising  spherical  cloud  can  be  described 
by 


H (t)  = R + Fd(R  (t)  - R ) 
sK  ' so  R s so' 


m 


(21) 


where  Fu  is  a constant  which  depends  upon  the  mass  entrainment  rate  a. 
K 

A reasonable  fit  to  F^  is 


!c  - 0 • 857 

5.31  a 

c -0.717 

5.31  a 


°-5  1 a 1 1 

1 < a < 2 


(22) 


It  is  a typical  effect  of  idealized  cloud  rise  models  that  the  ratio  of 
the  altitude  rise  to  the  radial  expansion  is  a constant.  Figure  9 
shows  an  example  of  the  agreement  between  the  analytic  rise  and  expansion 
fits.  Equations  19  to  21  , and  the  detailed  computer  calculations. 

We  next  formulate  analytical  expressions  for  the  rise  and  expansion 
of  the  munitions  clouds  by  the  appropriate  modification  of  the  idealized 
spherical  cloud  formulas.  The  initial  radii  of  the  dust  cloud  in  the 
track,  cross  track  and  vertical  directions  are  given  by  Equations  8 
and  11.  The  rise  and  expansion  of  the  dust  cloud  are  controlled  by  the 
state  of  the  atmosphere — primarily  wind  and  temperature  statistics. 
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The  Pasquill  Stability  Categories,  Reference  5,  are  commonly  used 
descriptors  of  the  atmospheric  state.  Another  parameter,  often  used 
in  similarity  equations,  is  the  atmospheric  turbulent  dissipation  rate, 
e (ergs  s g ”*).  The  Pasquill  parameters  arc  easier  to  measure  and 
arc  our  preferred  atmospheric  descriptors.  For  our  analytic  formulas, 
our  eventual  goal  is  to  correlate  a,  and  possibly  C^,  with  the  observed 
atmospheric  state.  In  the  related  program  of  analyzing  clouds  from 
countermeasure  smoke  devices,  the  cloud  behavior  is  being  correlated 
with  the  Pasquill  parameters. 

For  the  Dugway  munition  series,  tests  D1  and  D6  were  Pasquill 
Category  C (slightly  unstable);  the  other  tests  were  Category  D 
(neutral).  No  atmospheric  descriptors  were  available  in  the  pre- 
liminary Fort  Sill  data.  With  the  sparsity  of  data  at  this  time,  we 
will  not  be  able  to  correlate  the  analytic  parameters  with  the  atmospheric 
state.  Ne  expect  a fairly  high  air  entrainment  rate  due  to  the  initial 
and  induced  turbulence  during  the  rise  phase.  We  set  a equal  to  1 , a 
value  not  incompatible  with  our  limited  data  base.  We  set  C D equal  to 
0.8.  There  is  considerable  scatter  in  the  experimental  data;  but,  on 
the  average,  the  expansion  rate  in  the  horizontal  direction  is  about 
twice  the  expansion  rate  in  the  vertical  direction.  With  the  stated 
choices  of  a and  C^,  a reasonable  fit  to  the  experimental  data  is  to 
take  the  vertical  expansion  rate  equal  to  about  0.92  of  the  idealized 
spherical  cloud  radial  expansion  rate,  and  to  take  the  horizontal  ex- 
pansion rate  at  twice  the  vertical  rate.  Lacking  any  rise  experimental 
data,  we  simply  take  the  altitude  of  the  munition  cloud  eo’ial  to  the 
spherical  cloud  altitude.  Thus  the  munition  cloud  equations  during 
the  rise  phase  are 


ALTITUDE,  RADIUS  (it) 


where 


K./t)  = 


Rv(t)  = 


2KHt  + RIIO 


f2V  + RV0 


H(t)  = Ryjj  + 5.31(Rs(t)  - Rj) 


Rjj  = radius  in  horizontal  direction,  either  cross  track  or  along 
track 

Rv  = radius  in  the  vertical  direction  (m) 

Ky  = effective  vertical  turbulent  diffusion  coefficient 

= (0.92)2  Ks  = 0.715 v/fl V (m2  s"1)  (25) 

K..  = effective  horizontal  turbulent  diffusion  coefficient 


= 4 K, 


, 2 -1, 

(m  s ) 


Rj  is  given  by  Equation  8 , Rg  by  Equation  19  with  Rsq  = Rj.  R^q 

and  Ry0  are  calculated  from  Equation  11. 

We  are  assuming  that  the  effective  diffusion  in  the  horizontal  plane 
is  independent  of  direction.  Relative  to  the  center  o.f  mass  of  the  cloud 
the  diffusion  is  the  same  in  the  cross  wind  and  the  down  wind  directions. 
Studies  of  puffs  released  in  the  atmosphere  have  found  no  evidence  to 
contradict  this  assumption  (Reference  5 ) . Here  we  are  speaking  of 
the  hot  air  cloud  and  the  lightest  dust  particles  which  move  with  the  air 
velocity.  Heavier  dust  particles  will  lag  the  air  velocities  and  will 
fall  out,  leading  to  an  elongation  of  the  overall  dust  cloud  in  the 
wind  direction.  In  addition,  if  the  vertical  dimension  of  the  cloud 
is  large  enough,  the  shear  structure  of  the  wind  will  distort  the  cloud. 

For  comparing  the  analytic  formulations  with  the  experimental  data, 
we  use  the  vertical  and  cross  wind  test  dimensions.  In  the  down  wind 
direction,  the  length  of  the  primary  cloud  is  not  separated  in  the 
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reported  data  from  the  overall  length  which  includes  elongation  due  to 
fallout.  In  our  model  fallout  is  treated  by  considering  separate  volume 
elements  for  the  heavier  dust  particles.  There  is  no  doubt  which  is 
the  vertical  dimension  in  the  reported  data.  For  the  Dugway  data,  the 
width  (Y  dimension)  is  obviously  the  cross  wind  direction.  For  the 
very  preliminary  Fort  Sill  data,  it  is  not  readily  apparent  which  is 
the  wind  direction.  After  a survey  of  the  length  and  width  data  versus 
time,  we  conclude  that  the  width  dimensions  of  the  Tables  9 through  13 
are  cross  wind  dimensions  for  the  Fort  Sill  tests.  At  this  time  the 
cross  wind  dimension  cannot  be  determined  for  the  other  Fort  Sill  cases. 
For  the  test  comparisons  the  following  values  were  used  for  the  analytic 
calculations. 


Table  14.  Values  for  Dugway  analytic  calculations. 


Trial 

Yield  (lbs  TNT) 

Other  Parameters 

D1 

3.5 

f = 0.84 

^ = 1.18 

P 

(Dugway  elevation  « 4400  ft) 

Cj  = Cp  - Cy  - 1 

D2 

5.0 

D3 

6.75 

D4 

8.0 

D5 

15.0 

D6 

35.0 

I 


I 


. 


Table  15.  Values  for  Fort  Sill  analytic  calculations. 


Trial 


T-l 2 to  T-20 


Yield  (lbs  TNT) 


15.0 


6.75 


Other  Parameters 


f = 0.5 


— = 1.04 

P 


(Fort  sill  elevation  ss  1110  ft) 
CT  = 2 Cp  = 1 Cv  = 1/2 


The  yields  of  the  155  and  105  mm  Fort  Sill  rounds  were  not  specified  in 
the  preliminary  data.  Trials  D5  and  D3  of  the  Dugway  series  were  simu- 
lations of  155  and  105  mm  rounds,  respectively;  so  we  assume  the  same 
yields  for  the  Fort  Sill  rounds.  Figure  io  shows  the  Fort  Sill  height 
data  and  the  analytic  model  comparison  for  the  155  nun  round  test  T-3. 

The  fit  is  excellent;  in  fact  the  fit  is  better  than  could  be  expected 
from  the  statistical  variability.  There  are  no  width  data  for  the  155  mm 
round.  Figure  11  shows  the  height  comparisons  for  the  105  mm  Fort  Sill 
rounds.  This  figure  is  a good  illustration  of  the  variability  to  be 
expected  from  clouds  with  identical  yields.  The  mode)  is  again  a good  fit 
to  the  mean  of  the  test  data.  The  Fort  Sill  width  comparisons  for  the 
105  mm  round  are  shown  in  Figure  12  . The  T-16  cloud  is  considerably 
broader  than  the  model  prediction  at  later  times;  otherwise  the  model  is 
a reasonable  fit  to  the  mean  width  data.  Each  of  the  six  Dugway  tests 

has  a different  yield.  In  order  to  show  all  tests  on  the  same  figure, 

1/4 

we  have  normalized  the  dimensions  by  dividing  by  W . Except  for  the 
first  several  seconds,  this  normalization  will  reduce  all  model  predic- 
tions to  one  curve.  Figure  13  shows  the  Dugway  height  comparisons. 

There  is  such  a wide  variability  in  the  data  that  it  is  difficult  to 
assess  the  model  fit.  Figure  14  shows  the  Dugway  width  comparisons. 
Except  for  the  anomalous  fall  off  in  the  size  of  D1  after  20  seconds, 
the  model  is  a good  fit  to  the  mean  of  the  data. 
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Figure  11.  Comparisons  with  Fort  Sill  height  data  for  105-mm  rounds. 
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SECTION  4 

POST  RISE  PHASE 
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During  the  rise  phase,  turbulence  is  induced  in  the  dust  cloud  by 
the  differential  velocity  between  the  rising  cloud  and  the  ambient  air. 

As  long  as  the  induced  turbulence  is  the  dominant  diffusion  source,  the 
rise  model  equations  of  the  previous  section  are  expected  to  apply  to 
the  cloud  altitude  and  expansion.  As  the  cloud  rises  it  cools,  primar- 
ily from  mass  entrainment  but  also  from  adiabatic  expansion.  The  buoy- 
ant force  is  driven  by  the  temperature  differential  between  the  cloud 
and  the  ambient  air.  So  as  the  cloud  cools,  the  rise  velocity  decreases. 
Eventually,  the  cloud  stabilizes  in  altitude.  At  some  point  in  the  rise 
history,  the  induced  turbulence  becomes  so  weak  because  of  the  decreasing 
rise  velocity  that  the  ambient  atmosphere  turbulence  becomes  the  dominant 
diffusion  mechanism.  The  crossover  point  from  buoyant  driven  rise  to 
ambient  controlled  diffusion  can  occur  at  any  time,  depending  on  the 
ambient  atmosphere  conditions.  In  the  extreme  case  of  a very  high 
velocity  wind  with  strong  ambient  turbulence,  the  dust  cloud  may  be  torn 
apart  and  dissipated  near  ground  level.  In  this  extreme  case  the  cross- 
over time  occurs  essentially  instantaneously.  In  a calm  stable  atmos- 
spherc  the  crossover  point  probably  is  not  reached  until  cloud  stabili- 
zation at  late  time. 

The  model  equations  of  the  rise  phase  are  valid  for  the  dust  cloud 
rise  and  expansion  only  for  a limited  time.  We  first  consider  the  cloud 
rise  The  rise  model  equations  were  derived  assuming  an  idealized 
neutral  stability  temperature  profile.  For  comparisons,  the  computer 
rise  program  was  also  run  assuming  ideal  stable  and  unstable  temperature 
profiles.  The  following  is  a brief  description  of  these  atmospheric  models. 
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Assume  that  the  ideal  exponential  pressure  equation  obtains  for 
some  altitude  range  above  the  surface.  Then  in  this  altitude  regime 
H 


P = P c 
A o 


11 


SP 


dynes  cm 


(27) 


where 


-2, 


PA  = ambient  pressure  at  altitude  H (dynes  cm  ) 

_2 

Pq  = ambient  pressure  at  zero  altitude  (dynes  cm  ) 
Hgp  = ambient  pressure  scale  height  (m). 


The  adiabatic  lapse  rate,  I',  of  the  atmosphere  is  the  temperature  change 
with  altitude  due  to  adiabatic  expansion  (of  a rising  cloud  or  balloon). 
Using  Equation  16  we  can  write  the  adiabatic  lapse  rate  as 


r _ - t 

1 " dll  “ 


fill  J_ 

( Y PA  dH  j 


°K  m'1 


For  our  ideal  pressure  model 


(28) 


J_^A 

PA  dH 


H 


SP 


so  that 


f = - T 


(V  ) Hgp 


’K  m 


-1 


For  air  y = 1.4.  For  the  U.S.  Standard  Atmosphere 

Hgp  = 8.4  km  = 8400  in  at  sea  level 

T = 288  °K  at  sea  level 
o 

so  that  in  the  standard  case  at  sea  level 


-9.8  0K  km'1  = -9.8  x 10"3  0 


K in 


-1 


(29) 


(30) 


Let  the  temperature  of  a balloon  be  T , the  ambient  temperature  at  II  * 0. 
As  the  balloon  rises  and  expands  the  adiabatic  temperature  is 
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T - Toe 


(31) 


Y-l  II 

Y nSP  °K  . 


For  our 
to  be 


rise  model  atmosphere  we 


take  the  ambient  temperature  profile 
°K  (32) 


where 

= ambient  temperature  at  altitude  H (°K) 

To  = ambient  temperature  at  zero  altitude  (°K) 

Fg  = stability  factor. 

If  in  a given  altitude  range  the  ambient  temperature  profile  is  the 
same  as  the  adiabatic  temperature  profile,  then  by  definition  that 
portion  of  the  atmosphere  has  neutral  stability  with  respect  to  the 
temperature  profile.  If  the  ambient  temperature  falls  off  faster  with 
increasing  altitude  than  the  adiabatic  temperature,  the  atmosphere  is 
unstable.  In  a stable  atmosphere  the  adiabatic  temperature  falls  faster 
than  the  ambient  temperature.  For  our  model  atmosphere 


Fg  = 1 for  neutral  stability 
F&  > 1 for  an  unstable  atmosphere 
Fg  < 1 for  a stable  atmosphere. 

Let  a balloon  be  at  rest  at  altitude  H with  the  ambient  temperature. 
Let  the  balloon  be  displaced  slightly  in  the  upward  direction  and 
consider  the  buoyant  force  on  the  balloon.  In  a neutral  atmosphere  the 
balloon  temperature  drop  due  to  adiabatic  expansion  is  identical  to  the 
ambient  temperature  drop;  so  that  after  displacement  the  balloon  tem- 
perature is  still  equal  to  the  ambient  temperature.  The  buoyant  force 
on  the  balloon  after  displacement  is  zero.  The  balloon  is  "neutral"; 
after  displacement  it  has  no  tendency  to  move  either  in  the  up  or  down 
directions.  Next  consider  displacement  in  an  unstable  atmosphere. 

After  an  upward  displacement,  the  balloon  temperature  will  now  be  higher 
than  the  ambient  temperature.  There  will  be  an  upward  buoyant  force. 
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anil  the  balloon  will  continue  moving  in  the  upward  direction.  Once 
displaced  from  its  equilibrium  position,  the  balloon  is  "unstable." 

In  a stable  atmosphere  the  opposite  effect  occurs.  After  upward  dis- 
placement, the  balloon  temperature  is  lower  than  the  ambient  temperature. 
The  buoyant  force  is  negative  (downward)  and  the  balloon  is  returned  to 
its  previous  position;  the  balloon  is  in  stable  equilibrium.  Turbulence 
in  the  atmosphere  will  displace  parcels  of  air  from  their  original 
positions.  The  turbulent  induced  motion  can  be  amplified  in  an  unstable 
atmosphere,  suppressed  in  a stable  atmosphere,  or  unaffected  in  a neutral 
atmosphere . 

For  the  model  atmosphere  the  ambient  density  is  taken  as 


M P. 
o A 


PA  = R.T 


g cm 


-3 


o A 2.875  x 10  T, 


(33) 


What  are  the  characteristics  of  the  idealized  spherical  cloud  rise 
and  expansion  in  the  three  types  of  atmospheres?  Let  Pq  and  Tq  be  identi- 
cal for  all  three  atmospheres.  Then  the  initial  cloud  will  have  the 
same  size  and  temperature  in  each  atmosphere.  The  fractional  decrease 
in  cloud  temperature  due  to  adiabatic  expansion  is  the  same  in  each 
atmosphere  and  is  minor  compared  to  the  mass  entrainment  effect.  The 
temperature  profile  of  the  entrained  mass  is  the  primary  factor  control- 
ling the  cloud  rise  and  expansion.  The  early  rise  and  expansion  in 
each  atmosphere  is  essentially  identical;  it  is  only  after  the  cloud 
temperature  has  fallen  to  within  several  degrees  of  ambient  that  the 
differing  ambient  temperature  profiles  produce  differences  in  the  rise 
and  expansion  time  histories.  As  long  as  the  idealized  cloud  is  in  a 
neutral  or  unstable  atmosphere  it  will  continue  rising;  it  will  not 
stabilize  until  it  encounters  a stable  portion  of  the  atmosphere.  The 
cloud  rapidly  cools  to  within  a few  degrees  of  ambient,  then  asymptoti- 
cally approaches  the  ambient  temperature.  In  a stable  atmosphere,  the 
cloud  rises  and  cools  until  ambient  temperature  is  reached  and  thfe 
cloud  stabilizes.  Table  16  shows  the  stabilization  altitudes  for 

f* 

ideal  clouds  in  a stable  atmosphere  with  a temperature  stability  factor  of 
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Table  16.  Stabilization  altitudes  in  a stable 
atmosphere  with  F&  = -1. 


Hydro  Yield  (fW) 
lbs  TNT 

Stabilization  Altitude,  m 

a = 0. 5 

1 

2 

1 

90 

64 

46 

10 

146 

103 

77 

100 

239 

168 

130 

1000 

390 

278 

213 

Besides  using  the  model  atmospheres  described  here,  the  rise  program  can 
also  use  the  actual  measured  atmospheric  profiles,  if  they  are  known  for 
a given  munitions  test.  Then  the  stabilization  of  the  ideal  cloud  in 
the  test  site  atmosphere  can  be  computed. 

Thus,  given  sufficient  atmospheric  information,  we  can  calculate 
the  rise  and  expansion  characteristics  of  any  ideal  munition  cloud. 
Unfortunately,  the  behavior  of  real  munitions  clouds  can  differ  signifi- 
cantly from  an  ideal  cloud,  especially  at  latter  times.  Consider  an 
unstable  or  neutral  atmosphere  with  the  tropopause  (temperature  inver- 
sion) at  10  kilometers.  All  ideal  clouds  (and  any  real  adiabatic 
balloons)  will  rise,  however  slowly,  until  they  encounter  the  temperature 
inversion.  We  would  not  expect  a real  cloud,  say  from  a 1 lb  munition 
round, to  rise  10  kilometers.  In  a calm  windless  atmosphere,  the  1 lb 
cloud  might  retain  its  identity  long  enough  to  rise  several  hundred 
meters;  under  turbulent  windy  conditions,  the  cloud  might  not  reach  5 
meters  altitude. 

At  present,  there  is  little  experimental  data  for  guidance.  None 
of  the  rise  data  from  the  munitions  test  series  have  been  published. 
However,  the  related  problem  of  the  rise  of  a continuous  plume  from  a 
smokestack  has  been  studied,  both  experimentally  and  theoretically, 
for  many  years.  Cloud  and  plume  rise  in  the  real  atmosphere  is  a very 
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complex  process  and  it  is  not  surprising  that  no  complete  and  exact 
theory  has  emerged.  After  surveying  the  theoretical  and  experimental 
data,  Pasquill  (Reference  5 ) concludes  that  the  best  estimate  for  the 
stabilization  altitude  of  a continuous  plume  in  a neutral  or  near  neutral 
atmosphere  is  the  Briggs  formula  (Reference  6 ) 


pl/3  „*2/3 

AH  = 5 £; m (34) 

u 


where 


AH  = total  rise  from  altitude  of  emission  (m) 

F = stack  factor,  which  depends  upon  the  heat  content  and 

4 -3 

temperature  of  the  emitted  stack  gases  (in  s ) 


X*  = 0.45  F2/5  (B'^e-1)375  (m) 


B = ratio  of  Lagrangian  and  liulerian  time  scales  « 4. 

e = rate  of  dissipation  of  turbulent  kinetic  energy  per  unit 
mass  of  air  (joules  kg  * s *) 

u = average  wind  speed  (m  s *). 

Variations  of  as  much  as  50  percent  from  the  formula  are  observed  for 
real  plumes.  For  stabilization  in  an  unstable  atmosphere,  Briggs 
concludes  that,  pending  further  evidence,  the  neutral  formula  given 
above  is  the  current  best  estimate,  but  he  emphasizes  that  larger 
individual  fluctuations  must  be  expected.  In  terms  of  u and  e. 


H 


(u3  eV1/S 


(35) 


The  rise  and  stabilization  of  a discrete  munitions  cloud  are  expected 
to  differ  from  those  of  a continuous  plume.  We  take  the  stabilization 
altitude  in  neutral  or  unstable  conditions  to  be 


“max 


C„(fW) 


1/4 


(u3e2)1/S 


m 


(36) 
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where 


Cjj  = constant 

The  rise  rate  of  an  ideal  cloud  is  proportional  to  the  one-fourth 
power  of  the  hydro  yield.  The  u and  e dependencies  were  chosen  in 
analogy  with  the  plume  analysis.  The  constants  C..  can  be  set  if  suf- 
ficient  experimental  data  become  available  from  the  munitions  test 
series.  For  rise  in  a stable  atmosphere,  we  begin  with  the  stabiliza- 
tion altitude  predictions  of  the  ideal  cloud  rise  program.  We  then 
modify  these  estimates  to  account  for  the  average  wind  and  the  atmospheric 
turbulence.  Again,  final  forms  will  be  set  when  the  experimental  data 
become  available. 

The  Briggs  formula.  Equation  34,  and  our  analogy.  Equation  36,  are 
in  terms  of  the  atmospheric  dissipation  parameter  e.  This  parameter  is 
difficult  to  measure  and  is  probably  not  known  for  most  munitions  tests. 

A better  technique  will  be  to  correlate  the  stabilization  altitudes  with 
the  Pasquill  Stability  Categories,  when  sufficient  data  become  available. 

At  some  point  in  the  dust  cloud  rise  history,  the  ambient  turbulence 
becomes  the  dominant  diffusion  mechanism.  The  crossover  point  depends 
upon  the  ambient  atmosphere  parameters.  If  the  rising  cloud  encounters 
a temperature  inversion  and  stabilizes  at  a relatively  early  time,  then 
the  stabilization  time  is  probably  the  crossover  time.  For  other  clouds, 
the  crossover  time  will  occur  during  the  rise  phase. 

There  are  several  different  approaches  we  can  take  for  estimating 
the  atmospheric  diffusion  rates.  One  approach  is  to  use  similarity 
theory;  the  similarity  equations  utilize  the  dissipation  rate  parameter, 
e.  Another  method  is  to  correlate  the  diffusion  rates  with  the  Pasquill 
Stability  Categories.  We  will  develop  the  equations  for  both  the 
similarity  theory  and  the  Stability  Categories.  We  expect  to  adopt  the 
Stability  Category  method  in  our  final  model,  but  as  yet  we  do  not  have 
sufficient  data  to  make  a final  choice. 
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In  the  atmosphere  similarity  theory  and  experimental  observation 
indicates  that  the  cloud  diffusion  due  to  turbulence  has  the  following 
time  behavior  (References  7 and  8 ).  For  diffusion  in  the  horizontal 
direction 
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dt 
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For  diffusion  in  the  vertical  direction 


d oy(t) 


dt 


= 2 K_ 


(38) 


where 

2 2 

o = variance  of  the  cloud  concentration  (m  ) 

2 

o = initial  variance  of  the  cloud  in  the  horizontal  direction 
at  time  t (m^) 

tQ  = time  the  neutral  cloud  (i.e.,  cloud  of  same  temperature 
and  density  of  ambient  atmosphere)  was  introduced  into 
the  atmosphere  (s) 

e = rate  of  dissipation  of  turbulent  kinetic  energy  per  unit 
mass  of  air  (joules  g * s ) 

2 -1 

K = effective  constant  vertical  diffusion  coefficient  (ms) 

^ 2 

t,  = transition  time  from  t to  t behavior  for  horizontal  diffusion 
2/3  -1/3 


= C3  °H0 


(s). 
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and  C3  are  nondimens ions 1 constants  of  order  unity  resulting 


from  the  dimensional  analysis.  In  our  analysis  we  are  assuming  that  the 
relationship  between  the  cloud  variance  and  the  cloud  radius  is 


R(t)  = Cr0(t) 


(39) 


where  Cr  is  a constant,  generally  in  the  range  from  1 to  3.  Writing 
the  velocity  equations  in  terms  of  the  radius  we  have 
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where 


kva  ■ cl  h 
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The  subscript  A distinguishes  this  case  of  atmospheric  diffusion  of  a 
neutral  cloud  from  that  of  a buoyant  cloud. 


Let  the  crossover  time  from  buoyantly  controlled  diffusion  to 
atmospheric  diffusion  be  t£.  We  assume  the  radial  expansion  equations 
of  the  previous  section  apply  for  times  before  t£ , namely 
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(44a) 
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Define  the  cloud  radii  at  crossover  time 

= 2 Klltc  + RH0  m 

RVC  = 2 Vc  + RV0  m2  ‘ 


(45a) 

(45b) 


Then  ‘the  vertical  diffusion  after  crossover  time  is  taken  to  be  simply 


Ry (t)  = 2 KVA(t-tc)  + Rjc  m2  t > tc  . (46) 

The  vertical  dimension  of  the  cloud  is  continuous  over  the  crossover 
time,  but  the  expansion  velocity  changes  discontinuously. 
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For  t lie  horizontal  dimension,  we  assume  the  following  dust  cloud 
diffusion  behavior  after  crossover  time 
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t > t. 


arbitrarily  set  equal  to  unity.  The  fictitious  neutral  cloud 
parameters  R^qA  and  tQ  will  be  set  so  as  to  provide  size  continuity 
and,  if  possible,  velocity  continuity  at  crossover  time  tc>  The 
expansion  velocities  about  t£  are 


R^t) 

dt 


H 


0 < t < t 
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t^  is  the  time  the  atmospheric  diffusion  velocity  equals  the  buoyant 
rise  expansion  velocity.  The  crossover  time  is  taken  as 


tc  = MINttA»  *37]  s 


(50) 


where 


tgT  = stabilization  time  (s). 

When  the  stabilization  time  (tgj)  occurs  before  the  time  the  atmospheric 
diffusion  becomes  dominant  (tA),  the  velocity  is  discontinuous  at  the 
crossover  time  (t£) , with  the  velocity  lower  after  crossover.  This  case 
would  correspond  to  the  buoyant  cloud  encountering  a temperature  in- 
version and  stabilizing  fairly  early  in  time.  The  atmospheric  dif- 
fusion rate  is  still  relatively  weak  when  the  stabilization  ends  the 
buoyant  diffusion. 
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Thus  if  a neutral  cloud  of  half  the  size  of  the  munition  cloud  at  cross- 
over had  been  released  -t  seconds  before  detpnation  time,  both  the 
neutral  and  the  munition  clouds  would  have  the  same  size  and  expansion 
velocity  at  crossover  time  tc- 

For  t = t,.^  < t,  vie  choose  Rur.A  and  t so  as  to  maximize  the  ex- 
C ST  A nUA  O 

pansion  velocity  at  tc.  Thus 


^OA  = 


= *c  ‘ 2 


(A) 


Next  we  develop  the  equations  for  the  Stability  Category  diffusion. 
The  expansion  of  plumes  from  continuous  sources  (chimney  emission,  for 
example)  has  been  the  object  of  extensive  theoretical  and  experimental 
research.  The  research  indicates  that  the  diffusion  in  the  cross  wind 
and  vertical  dimensions  can  be  described  by 


o = a (Xn)  y 
y y v D' 


(55a) 
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(55b) 


where 

a ,a  ■ constants 

y * 

Xp  ■ downwind  distance  in  the  plume  at  which  the  o's  are 
evaluated  (m). 

For  the  individual  munitions  clouds,  time  instead  of  downwind  distance 
is  the  dependent  variable.  Since 


bl 


(50) 


where 


u = mean  wind  speed  (ms  ) 


we  take  the  standard  deviations  to  be 


Oy  = by(ut)  y 
°z  = bz^)  2 


(57a) 


(57b) 


Observations  indicate  that  the  instantaneous  (cloud)  an-J  continuous 
sources  have  the  same  power  law  dependence,  but  different  constant 
coefficients.  The  constants  and  power  exponents  are  functions  of 
the  Stability  Categories. 

For  matching  the  buoyant  controlled  rise  to  the  atmospheric  controlled 
diffusion  we  define  a fictitious  atmospheric  controlled  cloud  which  was 
introduced  into  the  atmosphere  at  time  tQ.  We  take  the  radii  of  the 
cloud  as 


RHA^  = Crby(u^t"tllO^ 


(58a) 


R„.  (t)  = C b [u (t-t  )] 
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(58b) 


The  t constants  will  be  defined  so  as  to  provide  radii  continuity  at 
the  crossover  time.  Again,  let  the  crossover  time  be  tc.  Radii  contin- 
uity requires 


tHO  = t 


t = t 
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and  R^  and  R^  are  given  in  Equation  45.  The  radius  of  the  munition 
cloud  is  given  by  Equation  44  before  crossover  time  and  by  Equation  5S 
after  crossover  time.  The  time  of  crossover  will  be  determined  from 


experiment;!]  observations  when  sufficient  data  arc  available. 


The  rise  phase  radial  expansion  behavior  is  a t ' variation. 

Similarity  theory  predicts  that  the  horizontal  radial  expansion  will 

1/2  3/2 

change  from  this  t variation  to  a t and  then  a t variation.  The 

Sy 

Stability  Category  predicts  at7  variation;  S^  varies  with  the  atmosphere 

stability  but  is  always  greater  than  one  half.  Examining  the  horizontal 

radial  expansion  data  in  Figures  12  and  14,  we  observe  no  noticeable 

1/2 

deviation  from  the  characteristic  rise  t behavior.  Either  the  cross- 
over time  has  not  been  reached,  or  the  data  reduction  technique  has  lost 
the  required  accuracy  to  show  the  deviation.  The  radius  is  taken  from 
photographic  data.  We  are  assuming  that  the  measured  radius  corresponds 
to  a constant  fraction  of  the  mass  penetrated  through  the  cloud  center. 

As  the  cloud  expands,  the  edges  become  wispy  and  faint  on  the  photographic 
image.  It  is  difficult  to  make  a consistent  measurement  as  the  cloud 
becomes  photographically  transparent.  Whatever  the  reason,  for  the 

present  we  have  no  munitions  test  data  for  the  expansion  behavior  after 
1/2 

the  t ' rise  controlled  phase.  Hopefully,  as  more  data  become  available 
from  the  munitions  tests  and  the  related  countermeasures  smoke  tests, 
the  behavior  during  the  atmospheric  controlled  phases  can  be  determined. 


SECTION  5 

DUST  MASS  LOADING 


When  a munition  round  detonates  at  or  below  the  ground  surface,  a 
quantity  of  soil  or  rock  mass  will  be  lofted  into  the  air.  Most  of 
the  lofted  mass  will  quickly  fall  back  to  earth  as  crater  ejecta  and 
fallback.  A smaller  fraction  of  the  mass  will  remain  aloft  in  the 
dust  cloud  for  long  periods.  Figure  15  shows  a sketch  of  an  idealized 
crater  formed  by  an  explosion.  The  total  mass  lofted  can  be  determined 
by  measuring  the  volume  of  the  true  crater. 


where 


total  mass  lofted  (g) 

average  density  of  soil  (g  cm 

3 

volume  of  true  crater  (cm  ). 


By  gathering  and  weighing  the  crater  fallback  and  ejecta,  the  mass  lofted 
into  the  dust  cloud  can  be  determined 


where 


dust  mass  lofted  into  the  cloud  which  remains  aloft  for 

long  periods  (g) 

mass  in  crater  fallback  (g) 

mass  in  crater  ejecta  (g). 


Terlccky  (References  9 and  10)  has  made  careful  mass  measurements 
on  5 crater’s  from  8 lb  TNT  detonations  and  1 crater  from  a 1000  lb  TNT 
detonation.  Not  surprisingly,  most  crater  data  consist  of  the  apparent 
crater  measurements  - radius,  depth  and  volume.  We  assume  that  the 
total  mass  lofted  and  the  dust  cloud  mass  arc  related  to  the  apparent 
crater  volume  by 


where 

3 

V.  = apparent  crater  volume  (cm  ) 

r\ 

= ratio  of  total  crater  mass  to  apparent  crater  mass 
F = ratio  of  cloud  mass  to  apparent  crater  mass 

w 

The  apparent  crater  volume  and  the  mass  ratios  depend  upon  the  explosive 
yield,  depth  of  burst  and  the  soil  type. 

We  first  consider  the  apparent  crater  volume.  The  volume  depends 
upon  the  soil  type.  We  intuitively  expect  craters  in  hard  rock  to  be 
smaller  than  craters  in  loose  soil.  We  expect  craters  in  dry  soil  to 
differ  from  craters  in  wet  soils.  The  soil  structure  will  affect  the 
crater.  For  a homogeneous  soil,  the  crater  volume  can  be  expected  to 
be  a smooth  function  of  the  yield  of  the  explosive.  In  a layered  soil 
this  may  not  be  true.  Consider  a soft  layer  overlying  a very  hard 
layer.  A large  yield  explosion  that  penetrates  into  the  hard  layer  may 
produce  a proportionately  smaller  crater  than  a small  yield  explosion 
whose  crater  is  entirely  in  the  soft  layer.  The  crater  volume  depends 
upon  burial  depth.  As  the  burial  depth  is  increased,  the  crater  volume 
increases  until  the  optimum  depth  is  reached.  As  the  depth  is  increased 
past  the  optimum,  the  crater  decreases  until  a great  enough  depth  is 
reached  so  that  the  overburden  completely  contains  the  explosive  gases 
and  no  ejecta  crater  is  formed. 

Vortman  in  two  review  articles.  References  11  and  12  , has  collected 


the  crater  measurement  data  from  a large  number  of  tests  of  near  surface, 
surface  and  buried  HE  detonations  in  soils  and  rocks.  The  depth  of 
burst  is  measured  from  the  charge  center  of  mass.  So  a half -buried 
spherical  charge  would  be  a surface  burst.  We  would  expect  from 
simple  similarity  considerations  that  the  crater  volume  would  scale  as 
the  first  power  of  the  yield.  For  large  HE  tests  (W  ^ 10^  lbs)  and 
nuclear  cratering  events,  the  scaling  does  appear  to  be  to  the  first 
power.  For  smaller  HE  craters  however,  the  scaling  power  appears  to 
be  slightly  larger  than  unity.  Figure  15Ashows  the  apparent  volume 
for  TNT  surface  bursts  for  soils  and  rocks.  A first  order  least  square 
fit  by  Vortman  to  the  test  results  for  craters  in  dry  lake  playa  at 


APPARENT  CRATER  VOLUME  (ftJ) 


Frenchman  Flat  at  the  Nevada  Test  Site  (NTS)  indicates  a scaling  of 

V = 0.21  W1*111  ft3 

A 

= 2.8  x loV'111  cm3  (63) 

The  Yucca  Flat  playa  also  appears  to  fit  this  scaling.  The  craters 
for  dry  clay  (Utah),  dry  sand  (Utah),  McCormic  Ranch  Playa  (New  Mexico), 
alluvium  (NTS)  and  limestone  (NTS)  appear  to  have  the  same  yield  power 
scaling  but  are  about  twice  the  volume  of  the  Yucca  Flat  playa  craters. 
The  Yakima  Firing  Center  basalt  craters  are  about  one-third  the  size 
of  the  Yucca  Flat  craters. 

For  buried  bursts  we  assume  the  same  power  scaling  (W*'***)  holds. 

We  use  similarity  theory  and  assume  that  the  depth  of  burst  normalizing 
parameter  is  the  yield  to  the  one-third  power.  Figure  16  shows  the 
apparent  crater  volume  data  and  a model  fit  for  TNT  bursts  in  alluvium 
at  area  10  (NTS).  Figure  17  shows  the  data  and  fit  for  alluvium  at 
Albuquerque.  The  Albuquerque  data  show  considerably  more  scatter  than 
the  NTS  data  and  seem  to  have  a narrower  distribution.  Figure  18 
shows  the  data  and  fit  for  dry  sand  in  Utah.  Generally  the  Utah  dry 
sand  is  very  similar  to  the  NTS  alluvium.  But  the  craters  for  the 
smallest  yields  (8  lbs  TNT)  are  about  three  times  larger  than  expected. 
This  departure  from  the  expected  yield  scaling  may  be  due  to  the  inhomo- 
geneity of  the  test  soil  with  depth.  Figures  19  through  24  show  the 
results  for  various  other  soils  and  rocks.  In  general  damp  soils  have 
larger  craters  than  dry  soils;  craters  in  hard  rock  are  smaller  than 
craters  in  soil. 

The  data  are  good  illustrations  of  the  uncertainties  associated  with 
cratering  and  dust  lofting.  Identical  yield  bursts  in  the  same  soil 
often  produce  craters  which  differ  by  a factor  of  2 or  more.  Different 
yield  bursts  at  one  site  may  or  may  not  scale  smoothly  depending  upon 
the  soil  characteristics.  Moreover,  live  munitions  may  crater  and 
scale  differently  from  these  static  test  events.  Until  data  become 
available  from  live  munitions  tests,  by  necessity  we  will  use  the 
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Figure  19.  Apparent  crater  volume  for  buried 
bursts  in  NTS  basalt. 
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Figure  20.  Apparent  crater  volume  for  buried 
bursts  in  Utah  sandstone. 
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Figure  21.  Apparent  crater  volume  for  buried 
bursts  In  Utah  dry  clay. 
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static  crater  data,  lor  predicting  crater  volumes,  at  a given  site, 
we  will  use  the  model  fit  for  the  soil  or  rock  type  that  is  closest  to 
the  type  of  soil  or  rock  at  the  given  site.  For  converting  volume  to 
mass  we  will  use  the  measured  soil/rock  density  at  the  site,  if  known. 
Otherwise  the  following  typical  average  values  will  be  used. 


At  present  we  have  little  direct  data  for  the  mass  ratios  and  Fc  of 
equations  62a  and  62b.  Table  17  shows  the  data  from  Reference  10. 

For  the  near  surface  and  surface  detonations,  of  events  1 and  5,  about 
25  percent  of  the  apparent  crater  mass  is  lofted  into  the  long  lasting 


Table  17  . Mass  fractions  for  McCormick  Ranch  playa. 


Event 

Yield 
(lb  TNT) 

Depth  of  Burial  (ft  pb 

Yield 

B9 

■n 
l| Hi 

1 

8 

-0.02* 

1.21 

0.27 

2 

8 

-0.04 

1.98 

3 

8 

-0.04 

-- 

— 

4 

8 

-0.04 

2.21 

0.044 

5 

8 

0 

1.35 

0.24 

6 

1000 

-0.04 

1.66 

0.099 

*A  minus  sign  indicates  above  ground  detonation. 


cloud.  The  fraction  falls  rapidly  with  altitude;  for  the  slightly  higher 
altitude  detonations  of  events  2,  4 and  6 the  cloud  mass  fraction  is 
only  about  4 percent  for  the  8 lb  charges  and  10  percent  for  the  1000 
lb  charge.  The  increased  cloud  loading  with  depth  of  burial  may  be  due 
to  additional  fracturing  and  soil  aggregate  breakup  by  the  explosive 
gases. 
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The  transmission  experiments  of  the  munitions  test  series  will 
provide  a source  of  indirect  data  for  dust  cloud  loading.  The  dust 
mass  penetrated  along  the  transmission  path  can  often  be  deduced  from 
the  test  data.  Utilizing  photographic  data  to  estimate  cloud  size, 
the  total  dust  mass  in  the  cloud  can  be  calculated.  When  sufficient 
munitions  test  data  become  available,  they  will  become  our  basic  data 
base.  Until  that  time  we  utilize  the  existing  static  charge  data. 


SECTION  6 

DUST  PARTICLE  SIZE  DISTRIBUTIONS 


For  reader  convenience  we  reproduce  the  first  five  pages  of  Section 
3 of  Reference  13  . The  equations  governing  the  two  most  common  and 
one  special  particle  size  distributions  are  discussed. 

The  two  most  common  probability  distributions  used  to  describe  par- 
ticulate (whether  dust,  smoke,  haze,  fog,  rain,  or  debris)  size  statistics 
are  the  power  law  and  the  log-normal  distributions.  The  power  law  prob- 
ability distribution  is 


(p  - 1)  a F 

:-(p-d  . a-(P-*) 

min  max 


a . < a < a 

mm  — — max 


where 


a = particle  diameter  (cm) 

Pp(a)  da  = fraction  of  the  particles  with  diameters  between 
n and  a ♦ da 

p = power  law  exponent 

a . = minimum  particle  diameter  in  the  distribution  (cm) 

min 

amax  * maximum  particle  diameter  in  the  distribution  (cm). 
The  log-normal  distribution  is 
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where 


= mean  particle  diameter  (cm) 

S ■ standard  deviation  parameter. 

The  two  distributions  given  are  probability  distributions,  ic. 
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Condensates,  such  as  water  droplets  or  the  particles  reformeJ  from 
the  weapon  vaporized  material,  are  generally  well  described  by  the  log- 
normal  size  distribution.  Some  typical  values  used  in  previous  parti- 
culate models  are  a = 1.4  * 10  5 cm  (0.14  pm),  S = 1.9  for  rccondcnscd 
ni  -4 

weapon  material  particulates  and  a = 5 * 10  cm  (5  lira),  S = 2.0  for 
water  droplets  within  a nuclear  cloud. 

Dust  and  crater  ejecta  particulates  are  generally  better  described 
by  a power  law  distribution,  at  least  for  particles  greater  than  a few 
hundredths  of  a centimeter  in  diameter.  Experimental  data  from  a number 
of  HE  (high  explosive)  and  nuclear  tests.  References  14  and  15  > indicate 
that  typical  values  for  the  power  law  exponent  are  p - 3.5  for  rock  and 
cohesive  soils  (such  as  clay  or  shale)  and  p - 4 for  loose  unconsoli- 
dated soils  (such  as  desert  alluvium  or  sand).  Fine  soils  may  have  p - 5. 

The  log-normal  distribution  is  well  behaved  mathematically  at  both 
limits  of  small  and  large  particle  diameters.  The  power  law  distribu- 
tion is  not  well  behaved  at  either  limit.  As  a . goes  to  zero  both  the 

min 

probability  distribution  P (amjn)  and  the  total  number  of  particles  .Vj.p 

go  to  infinity.  As  a goes  to  infinity,  the  probability  distribution 
niQ  x 


P(amaX)  goes  to  a well  behaved  zero,  but  NTp  goes  to  zero  for 


a unite 


Mr 


There  arc  several  techniques  employed  to  eliminate  mathematical  dif- 
ficulties at  the  endpoints.  The  most  common  technique,  and  the  one  we 
implicitly  assumed  when  we  defined  the  power  law,  is  simply  to  cut  the 

distribution  off  at  some  lower  and  upper  limits,  s . and  a . Another 

min  ntaX 

technique  is  to  assume  an  upper  limit,  a , hut  to  attach  a mathematic- 
ally well  behaved  tail  at  the  lower  end  of  the  distribution.  This  lower 
limit  tail  cpn  be  some  artificial  mathematical  expression  chosen  for 
computational  convenience,  or,  as  is  often  the  case,  can  be  a fit  to  the 
actual  experimental  data  in  the  small  particle  limit. 


DO 


SAI,  Reference  15,  reviewed  the  experimental  data  on  dust  particle 
size  distributions  from  nuclear  test  events.  They  then  analyzed  in  de- 
tail the  more  recent  (and  more  complete)  dust  size  data  from  IIH  tests. 
Their  conclusion  was  that  the  small  size  particles  have  a log-normal 
size  distribution,  while  the  larger  particles  have  a power  law  size 
distribution.  The  size  division  between  small  and  large  particles  oc- 
curs at  a particle  diameter  of  about  0.018  cm  (180  pm).  We  will  develop 
the  equations  for  the  hybrid  size  distribution  of  a log-noi;mal  small- 
limit  tail  attached  to  a power  law  distribution. 

The  hybrid  probability  distribution  is  given  by 


P„(a)  - 


ClPL(a) 


c2P_(a) 


0 < a < a . 

— — min 


a . < a < a 

min  — — max 


where  Cj  and  C 2 are  normalization  constants  which  ensure  that 


.(a)  da  = 1 


C.P.  (a  . ) = C,P  (a  . ,) 

1 nun'  2 pv  nun' 


ic,  the  constants  ensure  that  the  total  probability  distribution  is 
properly  normalized  to  unity  and  that  the  log-normal  and  the  power  law 
distributions  join  at  a jn.  and  C ^ are  given  by 


P (a  . ) 

C = P mln  c 
1 . I>  (a  . ) 2 

Lv  min' 


[ . /*"  TTN  V*.in>T‘ 

• rM^nvwJ  • 
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r 


r*  -it2 

'-“■7 tf 


dt 


(74) 


cumulative  probability  of  the  standardized  normal 
random  probability  distribution. 


For  the  hybrid  distribution, 

a 


. max  r nun  /•  max 

(a3),,  = f a3P(1(a)  da  = Cy  / a3l’L(a)  da  + C2  / a2P 
•'n  n 


a • 

nun 


max 


(a)  da 


a . 
min 


!s!»  \ 

1(*\  PuV-E^r"JtnS/  *c2<a  )i 


= C 


(75) 


and  the  total  number  of  particles  for  a given  mass  is 

“r 


N_..  = 


Z 


(76) 


For  the  hybrid  distribution,  the  fraction  of  the  total  number  of  particles 
in  the  power  law  segment  and  the  log-normal  segment  is 


fNP  = 

W|HJP  r 

nT„  = C2 

(77a) 

fNL  c 

°Wl  . c 

*ra  2 ' 

(77b) 

fractions 

of  the  total  dust  mass 

which  lie  in  the  two  segments  arc 

f ts 

*MP 

«r>p  - 

(78a) 
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fML  E 
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(78b) 

In  the  SAI  analysis  the  smallest  particles  measured  were  on  the  order 
of  about  0.002  cm  (20  pm).  The  smaller  particles  were  no  interest  for 
the  problem  being  analyzed.  However,  these  smaller  particles  may  be 
very  significant  with  respect  to  visible  and  infrared  transmission 
through  late-time  munition  dust  clouds.  We  require  improved  small 
particle  statistics. 

Patterson  and  Gillette,  Reference  16  , have  examined  the  size 
distributions  of  soil  derived  aerosols.  Under  conditions  of  heavy 
aerosol  loading,  they  conclude  that  there  are  two  distinct  size  modes 
present.  The  two  modes  are  log  normally  distributed.  The  smaller  mode 
has  a peak  near  1-ym  particle  diameter  and  seems  to  be  independent  of 
the  soil  type.  The  larger  mode  has  a peak  between  about  8 and  80  ym, 
depending  on  soil  type.  The  speculation  is  that  the  smaller  mode, 
labeled  mode  A by  Patterson  and  Gillette,  is  due  to  the  breaking  up  of 
the  loose  soil  aggregates  by  the  sandblasting  process  due  to  wind  erosion. 
The  larger  inode,  mode  B,  is  due  to  the  soil  aggregates  themselves.  In 
our  problem  of  dust  lofted  by  a cratering  explosion  we  would  expect  con- 
siderable fracturing  and  breaking  up  of  the  soil  aggregates. 

The  hybrid  distribution  measured  on  the  large  HE  and  nuclear  detona- 
tions appears  to  correspond  to  the  mode  B aerosol  distribution,  but  with 
a power  law  tail  for  the  explosively  ejected  larger  particles  (which  are 
not  present  in  the  wind  driven  aerosols).  The  mode  A particles  were 
not  measured  in  that  particular  series  of  tests,  but  have  been  measured 
in  some  late  time  nuclear  dust  clouds.  Dust  cloud  samplers  mounted  on 
aircraft  have  measured  the  size  statistics  of  dust  particles  in  nuclear 
dust  clouds  at  times  several  hours  after  detonation.  At  these  late 
times  particles  from  mode  A are  the  dominant  dust  constituent  in  the 
cloud. 

For  our  size  model,  we  adopt  the  two  mode  distribution.  We  assume 
there  are  two  distinct  size  groups  present,  the  mode  A and  mode  B 
particles.  We  take  the  mode  A particles  to  be  log  normally  distributed, 
and  the  mode  B particles  to  have  the  hybrid  distribution.  Table  18  gives 
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Table  18  . Size  probability  distribution  parameters. 


the  presently  recommended  values  for  the  probability  distribution 

parameters  (References  15  and  16  ).  Figure  25  shows  the  probability 

distributions.  The  mode  curves  shown  are  probability  curves;  that  is, 

* 

the  area  under  each  curve  is  unity.  The  fraction  of  the  lofted  dust  mass 
that  is  in  each  size  mode  can  be  determined  by  direct  size  measurements 
or  by  the  indirect  transmission  data.  For  wavelengths  of  a few  microns 
or  less,  the  two  modes  have  significantly  different  scattering  and  ab- 
sorption properties  per  gram  of  dust  mass.  If  we  have  an  estimate  of 
the  total  dust  mass  in  the  cloud,  then  the  mass  fraction  in  each  mode 
can  be  varied  until  a match  with  the  transmission  data  is  obtained. 


SECTION  7 

FALLOUT,  TRANSPORT,  AND  DIFFUSION 


The  particles  in  the  first  size  group  (see  Table  3 ) arc  small  enough 
to  have  negligible  fallout,  to  transport  at  the  wind  speed,  and  to  diffuse 
at  the  visible  cloud  rate.  Each  of  the  heavier  size  groups  will  have  a 
finite  fallout  rate,  lag  behind  in  the  wind  transport  and  diffuse  at  a 
slower  rate.  Let  the  wind  be  along  the  positive  (horizontal)  x axis  and 
take  the  z axis  positive  upward.  The  (horizontal)  y-axis  is  in  the  cross 
wind  direction.  Then  the  coordinates  of  the  centroid  of  the  1st  (smallest) 
dust  particle  size  group  are 


X 

o 

r+ 

►— * 

II 

m 

(79a) 

hTcWJi  = 0 

m 

(79b) 

[zc(t)]1  = H(t) 

m 

(79c) 

where 

Vw  = average  wind  speed  (m  s ) 

II (t ) = altitude  of  the  visible  cloud,  given  by  Equation  24  (m). 

In  some  cases  the  wind  shear  (variation  of  wind  velocity  with  altitude) 
will  be  so  large  that  an  average  wind  speed  over  the  rise  altitude  range 
has  little  meaning.  For  these  cases  we  use  the  more  general  expression 


[xc(t)]j  ./*  vw(t)dt  -/ 
0 0 


(H  (t)  w 


V (z) 

nr  v 7 


dz 


m 


(80) 


where 


l*c(t)]l  = 


dH(t) 

dt 


= rise  velocity  of  the  visible  cloud  at  time  t (ms  *), 


We  next  compute  the  z cooi'dinates  of  the  centroids  of  the  heavier 
particle  groups.  We  take  the  rise  (or  fall)  velocity  to  be  the  sum  of 
the  upward  buoyant  rise  velocity  and  the  downward  gravity  fall  velocity. 
The  upward  buoyant  rise  velocity  is  taken  to  be  the  same  for  each  size 
group  and  is  equal  to  the  visible  cloud  rise  velocity.  The  downward 
velocity  is  taken  as  the  terminal  velocity  under  gravity  fall.  Thus 
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^c(t)Ji  s [£c(t)1l  ' (Vi  m S_1 

(81) 

- (V^T  m 

(82) 

■ 

where 

i = 2 to  20,  the  index  of  the  particle  size  group 

• 

Izc(t)]^  = rise  velocity  of  the  centroid  of  size  group  i 

(m  s"1) 

(V  ).  = terminal  velocity  of  size  group  i under  gravity  fall 

z 1 _1 
(m  s ) 

[z  (t)].  = altitude  of  centroid  of  size  group  i (m) . 

C 1 

The  terminal  velocity  equation  of  a spherical  particle  falling  under 

gravity  in  still  air  is 

"*Pa  2 

4p  d"  V2  - S ’ 0 

P 

(83) 

where 

V 

= terminal  fall  velocity 

pa 

= ambient  air  density 

• 

PP 

= particle  density 

d 

= diameter  of  spherical  particle 

• 

g 

= acceleration  of  gravity. 

For  a 

spherical  dust  particle  we  can  use  (Reference  17  ) 

S - r * °-ta 

e 

(84) 

where 

is  the  Reynolds  number,  given  by 

paVd 

rc  - Hr- 

c pa 

i. 

(85) 

• 

and  u 

a 

is  the  air  viscosity.  Thus  the  velocity  equation  becomes 

I 

G7 
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V2  4 

3 v 


V - g = 0 


which  has  the  solution 


,//27u\‘  3Spnd  27pa 

V = V ± + 0 i_  (87) 

V \ Pad  / Pfl  Pad 

This  terminal  velocity  equation  is  more  familiar  in  the  limits  of  large 
and  small  dust  particles.  For  small  dust  particles  (d  < 0.01  cm) 

CD  * r and 

e 


which  is  Stokes  law.  lor  large  dust  particles  (d  > Oil  cm) 

Apjid 


V spacD 

which  is  the  dynamic  drag  equation.  We  assume  the  following  numerical 
values  for  air  near  sea  level 

p = 1.22S  x 10"3  g cm"3 

a -3 

Pp  *=  2.6  g cm 

g «=  980.6  cm  s 2 

p = 1.81  x 10-4  poises* 


With  these  values  the  terminal  velocity  equation  becomes 


/j'  3.9894^ 


, ,*6  , 3.9894 

+ 6.2438  x 1C  d ■ — 5 

d 


■1  (88) 


where  d is  in  cm. 


Now,  each  size  group  of  dust  particles  has  a range  of  particle  sizes 
in  it,  but  we  assume  the  whole  group  falls  with  one  characteristic  ve- 
locity. So  we  need  to  define  a characteristic  diameter  to  use  in  the 


terminal  velocity  equation.  We  make  the  simple  assumption  that  the 
characteristic  terminal  velocity  of  each  size  group  is  the  particle 
terminal  velocity  of  the  midpoint  diameter  of  the  size  group.  Table 
19  shows  some  terminal  velocities  in  sea  level  air.  The  buoyant  rise 


velocity  for  a cloud  with  a hydro  yield  (fW)  of  10  lbs  of  TNT  is  about 

n " 1 _ . . • r- ..  • 1 ^ 1 - 1 a.  f 


3ms  at  early  times.  For  this  case  particles  less  than  about  0.03  cm 
(300  pm)  in  diameter  would  be  carried  upward  for  some  distance  until 
the  buoyant  rise  velocity  slowed  sufficiently.  Larger  particles  would 
begin  falling  immediately. 


Table  19.  Terminal  velocities. 


Particle  Diameter  (cm) 

Terminal  Velocity 
(m  s-1) 

0 

0 

0.001 

0.0078 

0.003 

0.0702 

0.006 

0.276 

0.01 

0.718 

0.03 

3.20 

0.06 

5.49 

0.1 

7.51 

1.0 

24.9 

Next  we  consider  the  horizontal  transport  of  the  larger  particle 
groups  by  the  wind.  We  will  develop  the  equations  for  the  constant  wind 
case;  for  a varying  wind  the  velocity  equations  are  integrated  as  a func- 
tion of  time  (or  altitude).  Let 
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(89) 


lvD(t)li  = 


horizontal  difference  velocity  for  centroid  of 
particle  group  i 


[x  (t ) ] . = horizontal  velocity  of  centroid  of  particle 

c 1 -1 
group  i (m  s ) . 


Eliminating  the  gravity  term  from  Equation  86  , the  drag  force  per  unit 
mass  on  a spherical  particle  of  diameter  d is 


2 -1 
Fd  = CjVp  + C2VD  newtons  kg 


where 
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For  constant  V , 
w 


F„  = 


d[xc(t)J  dVD 


D dt 


dt 


Solving  for  VQ, 
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Integrating  with  respect  to  time, 


ii 
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j 
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[xc(t)].  = [xc(t)]1 


m 


For  the  characteristic  diameter  to  use  in  Equations  90  and  91  to 
evaluate  and  C ^ for  each  size  group,  we  again  take  the  midpoint 
diameter  of  the  range.  Table  20  gives  some  and  C ^ values  for  sea 
level  air  and  the  time  it  takes  for  the  particle  to  reach  90  percent  of 
the  wind  velocity  for  a 5 m s * wind.  As  would  be  expected,  the  small 
particles  rapidly  approach  the  wind  speed,  while  the  larger  particles 
take  progressive  longer  times  to  be  accelerated  by  the  wind. 


Table  20.  Horizontal  velocity  parameters. 


Particle  Diameter  (cm) 

C,  (in'.1) 

c2  (s’1) 

Time  to  Reach  90  Per- 
cent of  Wind  Speed  (s) 

0.001 

15.7 

1250 

0.00180 

0.003 

5.24 

139 

0.0155 

0.006 

2.62 

34.8 

0.0580 

0.01 

1.57 

12.5 

0.150 

0.03 

0.524 

1.39 

1.02 

0.06 

0.262 

0.348 

3.05 

0.1 

0.157 

0.125 

6.44 

0.3 

0.0524 

0.0139 

26.9 

0.6 

0.0262 

0.00348 

60.2 

1.0 

0.0157 

0.00125 

106.0 

The  larger  particles  will  expand  and  diffuse  at  slower  rates  than  the 
smaller  particles.  Following  Pasquill  (Reference  5)  we  have  approxi- 
mately 


wlicrc 


u = mean  air  velocity  (ms*) 

8 = ratio  of  Lagrangian  and  Eulerian  time  scales. 

For  8 we  take  an  average  value  of  4 (Reference  5 ).  In  the  model 

of  cloud  rise  and  expansion,  there  are  two  time  regimes.  In  the  first 
regime  buoyant  rise  and  expansion  dominate  atmospheric  diffusion;  in 
the  second  regime  atmospheric  diffusion  dominates.  In  the  first  time 
regime  we  take 


(K). 
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where 


K = effective  horizontal  or  vertical  diffusion  coefficients 


, 2 -1, 

(m  s J 


-1, 


Vjj  = initial  buoyant  rise  velocity  of  cloud  (ms  ). 


In  the  second  time  regime,  we  take 
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where 


2 R 

a = — =■  = variance  in  dust  mass  concentration  in  the  vertical 

Cr  2 

r or  horizontal  direction  (m  ) 
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Table  21.  Ratios  of  diffusion  parameters. 


Particle  Diameter 
(cm) 


0.001 

0.003 


0.006 


0.996 

0.938 

0.722 

0.228 

0.135 

0.099 

0.030 


Table  21  shows  the  ratio  of  the  effective  diffusion  coefficients 
or  the  variance  derivatives  for  or  equal  to  3 m s *.  Parti- 
cles larger  than  about  0.006  cm  (60  pm)  diffuse  at  the  full  cloud 
rate,  while  the  larger  particles  diffuse  at  decreasing  rates. 





SECTION  8 

SAMPLE  CALCULATION 


As  an  example,  we  will  compare  the  cloud  model  results  with  the 
experimentally  measured  data  for  Dugway  test  D4.  Figure  25  shows  the 
test  layout.  The  sight  paths  for  transmission  measurements  were  along 
rows  Q,  0,  and  M;  the  test  instruments  were  located  1.5  m above  ground 
level.  The  mean  wind  speed  was  4.7  m s * (measured  at  8 m)  blowing  from 
3'?  degrees  (about  6 degrees  off  vertical  with  respect  to  the  sight 
paths).  The  charge  was  8 lbs  of  TNT  buried  so  that  the  charge  top  was 
flush  with  the  ground  surface.  We  take  a wavelength  of  3.4  ym  for  the 
transmission  experiments. 

Figures  27  and  28  show  the  horizontal  and  vertical  radii  for  various 
particle  sizes.  Also  shown  on  the  figures  are  the  data  for  the  visible 
cloud  radii  which  we  are  assuming  correspond  to  the  zero  diameter 
particle  model  results.  The  agreement  is  quite  reasonable.  Figure  29 
shows  the  altitude  versus  time  for  various  particle  sizes.  We  have  no 
altitude  data  for  comparison.  Figure  30  shows  the  distance  the  various 
particle  sizes  are  carried  by  the  wind.  The  zero  diameter  particles 
travel  with  the  mean  wind  speed.  Note  that  we  are  ignoring  any  wind 
shear  and  simply  take  a constant  wind  speed  with  altitude.  Again  we 
have  no  data  for  comparison. 

Since  the  sight  paths  are  only  1.5  in  above  ground  level,  only  the 
larger  particles  (diameter  > ICO  ym)  will  fall  out  of  the  cloud  into  the 
sight  paths  in  the  time  it  takes  the  cloud  to  traverse  the  test  grid 
(~40  s).  Hence  we  can  ignore  the  minor  mass  in  the  smaller  mode  A 
particles  and  take  the  size  distribution  of  the  particles  to  be  the 
hybrid  distribution  for  unconsolidated  soils  of  Table  18.  Using  the 
cross  section  routines  given  in  the  appendix,  Table  22  gives  a printout 
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of  the  results  for  the  20  size  groups  of  Table  3.  The  index  of  refrac- 
tion chosen  is  a typical  value  for  soils  (References  21  and  22).  For 
particles  larger  than  about  30  pm  the  cross  sections  are  quite  insensi- 
tive to  the  index  of  refraction. 

The  charges  are  buried  so  that  the  top  of  the  charge  is  flush  with 
the  ground  level;  this  gives  a scaled  depth  of  burst  of 


D0B  - 0.04  m (lb  TNT)‘1/3 


W 


1/3 


(98) 


= 0.131  ft  (lb  TNT) 


-1/3 


If  we  assume  that  the  Dugway  soil  is  similar  to  NTS  area  10  alluvium  or 
Utah  dry  sand,  then  from  Figures  16  and  18 


W 


as  0.60  ft3  (lb  TNT)"1,111 
= 0.0170  m3  (lb  TNT)'1,111 


(99) 


_3 

Assuming  an  average  soil  density  of  1.5  g cm  , the  apparent  crater  mass 
is 


Ma  = 2.57  X 10  g 


We  have  a limited  set  of  data  for  the  ratio  of  the  total  mass  to  ap- 
parent mass.  The  average  value  for  the  5 cases  shown  in  Table  17  is 
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The  Gaussian  form  for  the  mass  distribution  is  particularly  con- 
venient for  horizontal  or  vertical  sight  path  integrations.  The  test 
grid  sight  paths  are  horizontal  and  arc  essentially  along  the  y-axis 
of  the  cloud.  Hence  the  mass  penetrated  (CL  value)  along  the  row  Q 
sight  path  for  the  volume  corresponding  to  one  of  the  size  groups  is 


f" 

MpQ  = P(X,  Y , Z)  cly 


e g m 

irrt  rr  ° 


2"Vy 


(101) 


where 

= X-ccordinate  of  centroid  of  particle  group  (m) 

a = standard  deviation  of  mass  concentration  in  the 
x 

X direction 

Xq  = X-coordinate  of  sightpath  (m) 

F = fraction  of  total  mass  in  the  size  group. 

M 

The  contribution  to  the  extinction  optical  depth  of  the  cloud  a ong  the 
sight  path  for  the  size  group  is 


x 


= MPQ  yE 


(102) 


where 

2 -1 

= mass  coefficient  for  extinction  (mg), 
b 

The  total  extinction  optical  depth  is  then  the  sum  of  the  contributions 
from  each  size  group 
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i=l 

and  the  transmission  is  given  by  the  usual  expression 


(103) 
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Wc  have  computed  the  total  optical  depth  for  rows  Q and  M.  The 
cloud  crosses  row  Q at  early  times  and  row  M at  later  times.  We  would 
expect  that  the  cloud  fallout  would  reduce  the  optical  depth  by  the  time 
the  cloud  reaches  row  M.  Figure  31  shows  the  model  results  and  the 
experimentally  measured  optical  depths.  Compared  to  the  test  data,  the 
model  results  are  lower  in  magnitude  and  narrower  in  time  distribution. 
The  general  behavior  of  both  the  model  and  the  data  is  what  we  would 
expect — the  early  line  of  sight  has  a higher  narrower  time  history 
than  the  later  line.  Considering  all  of  the  uncertainties  in  the  crater 
volumes,  mass  loadings,  fallout  rates,  mass  distributions,  particle  size 
statistics,  and  cross  sections  the  disagreement  in  the  magnitudes  is 
not  surprising.  In  fact  it  is  quite  encouraging  to  achieve  this  level 
of  magnitude  agreement  with  an  untuned  model  (the  model  parameters  have 
not  been  set  to  maximize  agreement  with  the  data).  The  disagreement  in 
time  spread  is  more  fundamental  however.  The  model  cloud  is  not  wide 
enough  near  ground  level.  The  expansion  rates  of  the  larger  particles 
may  be  depressed  excessively.  The  assumed  model  rise  rate  of  the  cloud 
may  be  too  large.  A slower  rising  cloud  would  both  increase  the  magni- 
tude of  the  model  results  and  increase  the  time  duration.  ’.Then  the 
test  rise  data  become  available,  this  effect  can  be  examined.  A more 
likely  explanation  is  that  both  the  extra  magnitude  and  time  duration 
are  due  to  a low-lying  dust  layer.  The  dust  model  implicitly  assumes 
that  when  a falling  dust  particle  impacts  the  ground  the  particle  stays 
there  (ie,  we  assume  a sticky  earth).  A large  dust  particle  impacting 
the  ground  can  break  up  or  scour  additional  dust  particles  from  the 
surface.  Smaller  dust  particles  impacting  the  ground  may  simply  bounce 
off.  The  surface  wind  can  carry  these  dust  particles  along.  The  test 
line  of  sig'ht  was  low  enough  (1.5  m)  to  be  in  the  surface  dust  layer. 
When  the  photography  data  become  available  we  should  be  able  to  tell 
which  effect  caused  the  lengthened  time  effect.  If  the  surface  dust 
layer  is  found  to  be  significant,  a dust  layer  model  can  be  added  to  the 
dust  cloud  model. 


Also,  tho  combustion  products  from  the  explosive  arc  not  yet  included 
in  the  dust  cloud  model.  Table  23  shows  the  yields  of  combustion 
products  per  kg  of  TNT  (Reference  23).  Carbon  has  a relatively  high 
extinction  coefficient  in  the  infrared.  The  modular  structure  of  the 
dust  cloud  computer  model  was  designed  to  include  any  number  of  additional 
particulates  besides  dust.  For  each  additional  particulate,  the  total 


1 


mass,  bulk  density,  size  distribution  and  complex  index  of  refraction 

. 

must  be  specified.  At  present  we  do  not  have  available  sufficient  data 
on  the  carbon  particulate  size  distribution.  When  these  data  are  obtained, 
. the  carbon  particulates  can  be  easily  added  to  the  dust  model. 


APPENDIX 

CROSS  SECTION  ROUTINES 


In  order  to  compare  model  predictions  with  the  experimental  trans- 
mission and  radiance  results,  we  need  to  calculate  the  scattering  and 
absorbing  properties  of  the  dust  particles  in  the  dust  cloud.  We  have 
developed  generalized  computer  routines  to  make  these  calculations.  In 
this  appendix  we  describe  the  program  routines. 

The  code  inputs  are: 

•.  Radiation  wavelength 

• Real  and  imaginary  parts  of  the  complex  index  of  refraction 
of  the  dust  particles 

• Density  of  the  dust  particles 

• Minimum  and  maximum  diameters  of  the  particle  size  group 
to  be  calculated 

• Type  of  size  probability  distribution  and  distribution 
parameters 

• Particle  shape  (whether  spherical  or  nonspherical) . 

The  code  outputs  are: 

• Average  scattering,  backscatter,  absorption,  and  extinction 
cross  sections  per  particle  in  the  size  group 

• Average  scattering,  backscatter,  absorption,  and  extinction 
mass  coefficients  per  gram  of  dust  in  the  size  group 

• Normalized  scattering  pattern  for  the  size  group 

• Number  of  particles  per  gram  of  mass  in  the  size  group 

• Number  fraction— the  ratio  of  the  number  of  particles  in 
the  size  group  to  the  total  number  of  particles  in  the  whole 
size  distribution 
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• Mass  fraction — the  ratio  of  the  mass  of  the  particles  in  the 
size  group  to  the  total  mass. 

Three  size  probability  distributions  are  built  into  the  routines — 
log  normal,  power  law,  and  hybrid.  These  three  distributions  are  dis- 
cussed in  Section  6 . Any  other  probability  distribution  including  a 
tabular  distribution  can  be  easily  added  to  the  program.  The  program 
is  set  up  to  calculate  the  outputs  for  a given  size  group  of  particles. 

This  size  group  can  be  a small  subset  of  the  whole  size  distribution, 
as  in  the  case  of  one  of  our  20  dust  particle  size  groups.  Or  the 
size  group  can  include  the  whole  probability  distribution. 

The  standard  Mie  calculation  routine  for  single  spherical  particles 
is  a new  routine  which  was  recently  developed  for  the  WOE  code  (Reference 
18  ).  The  routine  uses  a stable  backwards  recursion  technique,  utilizing 

Lentz's  (Reference  19  ) continued  fraction  algorithm  to  evaluate  the 
beginning  term.  The  backward  recursion  technique  does  not  suffer  from 
the  instabilities  and  inaccuracies  of  the  forward  recursion  techniques 
of  many  earlier  Mie  codes.  The  present  Mie  routine  is  compact,  fast 
running,  and  accurate  for  any  combination  of  wavelengths  and  particle 
sizes. 

Dust  particles  are,  in  general,  irregular  and  nonsphcrical . The 
spherical  Mie  solutions  are  normally  used  to  calculate  the  scattering 
and  absorption  properties  of  any  nonspherical  particle  that  is  not 
planar  (ie,  flat,  such  as  a snowflake).  The  rationale  is  that  the  non- 
spherical particles  are  randomly  oriented  so  that  on  the  average  their 
scattering  and  absorbing  properties  are  the  same  as  spherical  particles 
with  the  same  volume  or  cross  sectional  area.  For  total  scattering  and 
absorbing  cross  sections,  the  results  are  generally  quite  accurate  for 


the  nonspherical  particles.  However,  the  details  of  the  scattering 
pattern  are  different  for  the  spherical  and  nonspherical  particles. 

The  back  and  side  scatter  from  a distribution  of  nonspherical  dust 
particles  can  differ  significantly  from  that  of  an  equivalent  spherical 
distribution.  There  are  various  techniques  available  for  the  calcula- 
tion of  the  scattering  pattern  of  nonspherical  particles.  Pinnick,  et  al, 
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in  Reference  20  have  shown  that  a simple  modification  to  the  standard 
Mie  solutions  yields  scattering  pattern  results  for  nonspherical 
particles  that  are  in  good  agreement  with  experimental  data.  The 
modification  eliminates  the  effects  of  the  surface  waves  that  are 
present  in  scattering  by  spherical  particles  but  are  absent  in  scat- 
tering by  nonspherical  particles. 

A preliminary  version  for  nonspherical  Mie  calculation  based  on  the 
Pinnick  method  has  been  incorporated  in  the  routines.  The  outputs  from 
the  nonspherical  calculations  are  in  general  agreement  with  the  examples 
given  in. Reference  20.  But  further  development  is  needed  before  the 
nonspherical  option  can  be  used  with  confidence  in  all  cases.  The 
routines  described  here  were  developed  as  a tool  to  use  for  comparing 
measured  transmissions  with  the  calculated  transmissions  through  the 
model  dust  cloud.  The  development  of  the  dust  cloud  model  is  the 
primary  purpose  of  the  present  effort.  Thus,  minimum  effort  was  devoted 
to  the  Mie  routine  development.  The  nonspherical  option  can  be  improved 
in  some  future  effort. 

A new  variable  step  length  algorithm  has  been  developed  for  the 
integration  over  particle  sizes.  The  algorithm  accounts  for  the  rate 
of  change  of  integrand  magnitude  with  step  length  and  the  contribution 
of  the  step  to  the  total  integral.  A rapidly  changing' integrand  requires 
smaller  steps,  while  a portion  of  the  integral  that  contributes  little 
to  the  total  can  be  calculated  with  larger  steps.  The  algorithm  is 
designed  to  produce  maximum  accuracy  with  minimum  computational  time. 

The  equations  solved  in  the  routines  are  the  following  The  mean 
cross  section  per  particle  is 
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- ,/al  1“  Q(a)P(a)da  , 
a = r cm  per  particle 
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where 


1 


* average  cross  section  per  particle  for  scattering, 

2 

backscattcr,  absorption,  or  extinction  (cm  ) 

■ particle  diameter  (cm) 

■ minimum  diameter  of  particle  size  group  (cm) 

a 2 * maximum  diameter  of  particle  size  group  (cm) 

Q(a)  ° Mie  efficiency  for  scattering,  backscatter,  absoiption, 
or  extinction  for  a particle  of  diameter  a 
P(a)  = size  probability  distribution  (log  normal,  power  law, 
or  hybrid,  see  Section  6 for  equations). 

The  normalized  scattering  pattern  is 


SN(9)  = 


°SCA  •'a. 


2 TTO  2 

QSCA(a)P(a)S(6,a)da 


(A- 2) 


where 


OgcA  = average  cross  section  per  particle  for  scattering, 
given  in  Equation  A-l  (cm^) 

Qsca(s)  = Mie  efficiency  for  scattering 

S(0,a)  = normalized  scattering  pattern  at  scattering  angle  0 

for  a single  particle  of  diameter  a 
= normalized  scattering  pattern  for  the  size  group. 

The  scattering  patterns  are  normalized  so  that 
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S(0)d«  = 1 


(A-3) 


Hence  an  isotropic  scattering  pattern  is  identically  unity  for  all 
angles.  The  coefficients  per  gram  of  material  in  the  size  group  are 


Nt  O 


2 -1 
cm  g 


(A-4) 


where 
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p = mass  coefficient  for  scattering,  backscattcr,  absorption, 
or  extinction  (cm^  g *) 

Nt  = total  number  of  particles  per  gram  of  material  in  the  size 


group  (g"1), 


The  integration  strategy  is  the  following 

,'a. 


In 
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j=i  4. 
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where 

1(a)  = integrand  of  the  integral  of  interest. 


The  interval  from  to  a2  is  broken  up  into  a number  of  subintervals. 
For  each  subinterval  we  assume  that  the  integrand  can  be  described  by 
a power  law 
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where 
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Then  the  contribution  to  the  total  integral  of  this  subinterval  is 

I (a)da 
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The  contributions  of  the  subintervals  arc  calculated  and  summed  until 
cither 

(1)  the  diameter  limits  of  the  size  group  interval  are  reached, 

(2)  or  the  integral  converges  to  its  final  value  (the  remaining 
unsummed  subintervals  are  negligible). 


The  step  size  is  chosen  as  the  e-folding  distance  for  the  integrand 

1 


(Aa) 

With  this  step  size 
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where  the  plus  sign  is  for  positive  x and  the  negative  sign  for  negative 
x.  In  addition,  for  those  cases  where  the  integrand  is  decreasing  so 
that  the  subinterval  contributions  are  decreasing,  the  step  size  is  in- 
creased. The  increase  is  taken  as  the  ratio  of  the  previous  contribution 
to  the  last  contribution 
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